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Chapter 1 
Introduction 



The binary stars are crucial for our knowledge about the universe. Especially eclipsing 
binaries provide us an unique insight to the basic physical parameters of the stars, stellar 
clusters, interstellar medium and galaxies. They are excellent distance indicators. We are 
able to learn more about the matter composition of the stars, about their evolution status, 
or the presence of planets or other components in these systems. 

The very first task is the data acquisition, because only with precise input data is one 
able to get precise results. In last few decades mainly due to excellent satellite observatories 
(and not only in the visible part of the spectrum) our knowledge of them has rapidly grown. 

Regarding the astrometry, there is still decreasing the number of observations of the 
wide pairs. On the other hand, due to the new interferometers, which could resolve 
the milli- and micro- arcsecond angular distances, the observable semimajor axes of the 
astrometric binaries are still decreasing. Unfortunately, most of the systems analyzed 
below have the angular size of the astrometric orbit from 1 arcsec down to 100 mas, which 
is beyond the limits for the modern multi-aperture interferometers. And the lack of recent 
observations lead to the low accuracy of the results. 

Another approach is photometry and the classical observation of minimum light. Due 
to the large "baseline" of observers in our country and the interest of amateur astronomers, 
the number of these observations is growing very rapidly and the cooperation between pro- 
fessional and amateur astronomers is very intensive. Many of the observations of minimum 
timings used in this study came from amateur astronomers and these measurements are 
as accurate as from the professional observatories. Thanks to the large minimum times 
data set we are able to analyze many of the eclipsing binary systems for their long-term 
period variations. 

The whole thesis is divided into several parts. In the first one is presented the theory 
needed for the analysis of multiple stellar systems by photometric and astrometric tech- 
niques, description of such systems and some limitations which have to be considered. In 
the second part are introduced several systems which show apparent period changes in 
their O — C diagrams. And in the third one are the systems analyzed by photometry and 
astrometry simultaneously. Also the catalogue of other suggested systems for simultaneous 
analysis is included in this chapter. This is the crucial part of this thesis. The method 
itself is introduced in the chapter [2] and the results are in the chapter HI 
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Chapter 2 
Theory 



2.1 Binaries 

Most stars are found to be members of binary or multiple stella r systems. A recent 
analysis of a large set of close binaries f Pribulla fc Rucinskil 2006 ) indicates that even 



most of the binaries are in multiple systems. According to this study about 59 % of 
northern-hemisphere contact binaries are members of multiple stellar systems. The sample 
of binaries was analyzed very precisely, some distance-independent techniques were used, 
and the selection effects were also discussed. 

During the last decades, many of the observational techniques have become so effec- 
tive and precise that it is possible to discover low-massive stars, brown dwarfs, or even 
exoplanets in eclipsing binary systems (hereafter EB). 

The astrometric binaries are a special subset of binaries, where the individual com- 
ponents could be resolved into separate stars, which means the angular separation of the 
components has to be above a certain limit (a function of the telescope aperture and 
technique used). Because the individual components in the system were discovered at 
different time epochs, they were marked by different labels. Most common is the use of 
A-B-C. . . sequence, which means that the component "A" was discovered as the first one, 
then was discovered the second one " B" , and after then " C" , etc. Sometimes one compo- 
nent was resolved to be a double, so it turns A-B Aab-B. Mostly A component is the 
brightest one. Here comes the problem with the hierarchy of such a system. Sometimes it 
is so-called hierarchical-type Aabc-B system (where the B component is far away from the 
triple Aabc, where the c component could be far away from the ab double), and sometimes 
is is so-called trapezium-type Aab-Bab system ( two doubles Aab an d Bab far away from 
each other), for the detailed description see e.g. iDocobo &: Andradd fl2006h . 



In the present analysis another approach was used. It is based on the physical properties 
and gravitational bounding of the stars. The numbers 1 and 2 were used for the primary 
and secondary component of the eclipsing pair (these are not spatially resolvable) and 
index 3 as a label for the third, distant, component which is astrometrically observable. 
Sometimes an additional component, the fourth, is suggested. In this study only the 
hierarchical systems are analyzed. It means that the higher the number of the component, 
the bigger the semimajor axis (the distance of the component is larger than the previous 
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one) ai <^ a2 <^ 03 ^ . . .. Because here we deal only with the relative astrometry (position 
angles and angular separation) it should be more precise to write au <^ ^12-3 ^ • ■ ■• 



2.2 Principal methods for analyzing the EBs 

Despite the fact that the methods introduced below are the most important and useful 
ones, only a short description was presented here. These methods are not the essential ones 
for this thesis. Two other methods (the astrometry and the analysis of times-of-minima 
observations) are more crucial and are described in detail in the next section. 



2.2.1 Spectroscopy 



Observing the spectra of stars is perhaps the most time-consuming activity performed 
at astronomical observatories all around the world. However, obtaining the spectrum 
of the star is also the most powerful tool for deriving the relevant parameters of the 
star. Thanks to different techniques of spectral analysis (spectrophotometry, line-profile 
analysis, disentangling, etc.) one can model the stellar atmosphere, derive the orbital 
parameters, or discover the surf ace structures. For a brief introduction to the topic and 
overview of the methods see e.g. iHilditchI (120011 ). 

Besides spectral classification, the modelling of radial velocity curves (hereafter RV) 
is one of the oldest methods of spectral analysis. If only one component of the binary 
is observable, we deal with a so-called SBl-type, and if both components are evident, it 
is called SB2-type binary. From both RV curves, we could derive many parameters of 
the relative orb i ts of the stars in the system (see Table ETT] for the parameters and e.g. 
Wilson fc Sofial (1l976l ) for the na e thod) . Cross-correlation methods are now often used for 
analyzing RVs, see e.g. iMcLeanI fll98lh. 



Spectral disentangling (see e.g. 



Hadravdll995l ) of composite spectra into the separate 
ones helps us to model the particular star in the system. Theoretical stellar models could 
be compared to the observed ones and one could study the physical conditions in the 
atmosphere of the star, temperatures, pressures, rotation, stellar wind or the chemical 
composition. 



Another quite new technique is Doppler profile mapping (see e.g. iRice et al.lll989l ). 
With this technique one is able to discover the dark and bright (cool and hot) spots on 
the surface of the st ar, as well as its rota tion, or the evolution of the spots. Doppler 
tomography (see e.g. iMarsh fc Horndll988l ) is able to detect similar structures and effects 
in accretion discs in interacting binaries. Evidence for discs, jets, or outflows could be 
observed in the precisely measured spectrum of the star. 



Table 2.1: The scheme of directly derivable entities. 
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Some comments: The table shows which entities could be derived when one has measurements of a particul ar type. The table 
i s div ided into two parts. In the first one (the left top corner) are the parameters and methods adopted from iKallrath &: Milone 
(|l999l b page 142. In the second part are the parameters of the third-body orbit which could be derived from the methods used in 
this thesis. 'LC stands for the light curve, 'RV for the radial velocity curve, 'SBl' and 'SB2' for the types of the spectroscopic 
binaries, 'LITE' for the light-time effect, 'Apsid.' for the apsidal motion and 'Astrom.' for the astrometric orbit, respectively. 
In the parameters column ai and 02 denote the semimajor axis of the primary and secondary component in the EB, i for the 
inclination of the EB orbit, Mi and M2 masses of the primary and secondary, Ri and R2 radii of components, Ci and £2 for the 
bolometric luminosities, d for the separation between the components, e for the eccentricity of the EB orbit, u for the argument 
of periastron of the orbit, 7 is the systemic velocity of the EB pair, q is the photometric or the spectroscopic mass ratio, Li and 
L2 monochromatic luminosities in a specified passband, gi and §2 gravity darkening coefficients, Ai and A2 albedo coefficients, Fi 
and F2 rotation parameters, Fi = xi and X2 limb-darkening coefficie nts, ^3 the third ligh t and T2 temperature of the secondary 
(Ti fixed), respectively. Only one parameter was added compared with IKallrath Sz Milond . the apsidal motion rate uj = and 
parentheses mean that the parameter is computable only if some apsidal motion is presented. Sometimes it is difficult to compute 
the parameters, or the ability to compute them depends on the type of the EB. For example the photometric mass ratio could 
be computed only when the EB is contact or over-contact type (that is the reason why there are the parenthesis). The question 
mark indicates that the temperatures could be derived from the spectrum of the star. The parameters from the second part are 
explained in the text, D stands for the distance to the system. 
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Sometimes also an additional component is observable in the spectra. The spectral 
lines of the third component typically remain at a fixed wavelength, while the lines of the 
binary components move in agreement with the actual orbital phase of the EB. In the 
more favourable case the third lines also move very slowly according to the third-body 
orbit. In some others only the barycenter of the whole system is moving very slowly in 
agreement with the phase of the third-body orbit (see Eq. 12.71 in chapter l2.3.4p . Such 
effec ts are hardly observable and their final detection is questionable in most of the cases 



see 



Mayerl (120041 ) for the list of such systems and the subsections 14.21 and 14.41 below) . 



2.2.2 Photometry 



Modern photoelectric and CCD photometry is a very powerful tool for studying the eclips- 
ing binaries. Thanks to a very wide network of amateur astronomers with their CCD 
cameras all around the world, it is nowadays quite easy to observe the whole light curve 
of a particular eclipsing binary in the standard filters and perform a detailed analysis of 
the system. 

From the light curve (hereafter LC) one could derive many useful parameters of the 
eclipsing binary itself. The comparison between the parameters derivable from the LC 
solution, from the RV fitting and from the techniques described below is shown in Table 

o 

From this set of parameters, together with the radial velocity curves, one could obtain 
a complete set of parameters describing the orbit of the individual components in the 
EB, as well as the basic physical parameters of both components in absolute units. The 
masses, the radii, the luminosities and the temperatures could be calculated, if the precise 
photometry together with both radial velocities w ere carried out (see Table [2TT1) . The limb- 
darkening coefficients (see e.g. Ivan Hammelll993h . as well as gravity brightening (see e.g. 
Lucylll967l ) and albedos (see e.g. lR,ucihskilll969l i of the components could be also calculated 



from the model. The most commo n codes for the EB modelling are the Wilson-D evinney 
flWilson fc Devinnevl[l97lh . FOTEL flHadravall2o"o3 l. Linnell's model (lLinnelllll98i ). LIGHT 
EilL,1979. ). etc. 

Modern advanced tools for analyzing the light curves of EBs are very powerful for 
discovering the surface structures and their evolutio n. Dark or hot sp ots on the star 
surface could be implemented into the model (see e.g. |Poe fc Eatonlll985l ). 

Precise photometry of the binary in different filters could also reveal the third light 
from a distant component in the LC solution. The third light I3 could be the indicator or 
the proof for the presence of a third body suggested by some independent method. 



2.2.3 Additional techniques 

Besides photometry and spectroscopy other methods also play a role in our modern knowl- 
edge about EBs. The invisible part of the spectrum is also used to study these objects. 
Some EBs have been successfully identified to be radio- or X-ray active. Such systems 
are mostly the complicated ones, currently in an unstable evolutionary status, undergoing 
rapid mass transfer, having an accretion disc or be a member of a cataclysmic variable. 



J. ±±j^Ky ± L J. 



Another appro ach is the use of polarimetry (see e.g. lHal]|ll949l ). or magnetometry (see 



e.g. iMarcyl Il984l ) to study EBs. These two modern techniques are also very powerful 
and could reveal some properties of these systems which are otherwise undetectable, for 
example a disc or stream in the system, large coronae of one of the components, or magnetic 
fields. 



2.3 O — C diagram analysis 

The periodic behavior of eclipsing binaries can also be studied. Long-term changes of 
the binary period could be caused by various effects and these effects could be analyzed 
thanks to the large database of times-of-minima observations which cover more than a 
century in many cases. Especially due to amateur astronomers with their CCD cameras 
in the recent few years the number of times of minima is growing rapidly, because these 
observations are much easier to obtain than the photometric observation of the whole LC. 
The accuracy of such observations are sometimes not very good (mainly the old ones), so 
the an alysis is also problematic. The topic was discussed in detail for example in Isterkeni 
fl2005h . 

The O — C diagram in our case is a special plot, where the x-axis is either the epoch 
(number of cycles relative to JDq) or the Julian date,, while the y-axis is the difference 
Observed minus Calculated. Here in the case of the eclipsing binaries this means the 
difference in times of minima in the particular system, expected minus predicted moment 
of minimum light. 

The linear ephemeris of the binary indicates the minimum time JD after E cycles {E 
is the epoch number) since the initial time of minimum JDq occured. One could write 

JD = JDo + P ■ E, 

where P is the period of the eclipsing binary. The value O — C is therefore defined as a 
difference 

O - C = JD - JDo - P ■ E. 

Finding new revised linear ephemeris of the binary could be also done with these 
equations. One has to minimize the sum 

N 

^(O-C).^-O 

i=l 

with respect to the parameters {JDq, P) over the whole parameter space and is the 
number of times of minima used. 



2.3.1 Constant period 

One would expect a constant period and therefore a linear trend in the O — C diagrams 
in most of the eclipsing binaries. If the assumed period of the EB is lower than the right 
one, the linear trend in the diagram is increasing, while if the proposed period is higher 
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than the right one, the times of minima are decreasing in the diagram. A collection of 
many O — C diagra ms of EBs is for example in Kreiner's Atlas of O — C diagrams, see 



Kreiner et all ((20011). 

Sometimes O — C diagrams show changes in period which seems to be abrupt. Such 
"jumps" in the diagram could be caused by sudden period changes. These are often 
explained as a mass ejection and/or transfer from one star to another. Sometimes this 
explanation is used only due to poor coverage of the abrupt change in the diagram, some- 
times is confirmed by another independent method. 



2.3.2 Mass transfer 

Another phenomena which could be studied in the — C diagram is the parabolic behavior 
of times of minima. If the data set is large enough, in some binaries it is possible to identify 
a parabolic (increasing or decreasing) trend, which means that the period of the binary is 
steady increasing or decreasing. In O — C diagram it means 

JD = JDq + P ■ E + q- E"^, 

where q is the quadratic term coefficient. It is possible to show that this additional term 
could be caused by the mass transfer within the binary. 

There are two basic kinds of mass transfer between the components. The first (and 
the simplest) one is conservative mass transfer. In this case the total mass of the binary 
as well as its total orbital angular momentum are conserved in the system. The other one 
is nonconservative mass transfer, where these quantities are not constant for the whole 
system. A few mechanisms could cause this latter case: stellar wind, Roche-lobe overfiow, 
or a sudden catastrophic event. This kind of mass transfer is probably more often in the 
nature. 

Concerning the former case, from the quadratic term coefficient q also the conservative 
mass transfer rate could be derived 

M1M2 P 2qMiM2 



3(Mi - M2) P 3P2(Mi - M2) ■ 
A brief i ntroducti o n to t he topic and derivation of the equations for both cases could be 



found in iHilditchI (120011 ) . page 162. The typical values of mass transfer rate are circa 



10"^ - 10-9 Mo/yr. 



2.3.3 Apsidal motion 

The next effect which could be studied only on the basis of the — C diagram analysis 
is the apsidal motion. In some eccentric binaries the line of apsides of the orbit of such a 
system is moving in space. One has to take into consideration two different periods, the 
sidereal one P, and the anomalistic one Pa, which are in relation 

Ps = Pail-i:o/2n). 



The quantity uj is the apsidal motion rate. 



Such an effect is easily detectable in the O — C diagram of the binary. Both primary 
and secondary minima are being periodically shifted from the linear ephemeris and also 
against the other one (the primaries and secondaries are in anti- correlation). The necessary 



equati ons for the apsidal 0—C diagram analysis are presented in lGimenez fc Garcia-Pelayo 



(119831). 



2.3.4 Light - time effect 



Large set of times of minima of the EBs could be also used for discovering the additional 
component (s) in these systems. With the light-time effect (hereafter LITE) analysis one 
can suppose a presence of another component (s) in the system only by analyzing the 
times of minima and their long-time behaviour. The motion of the EE around the com- 
mon barycenter causes apparent changes of the observed binary period with a period 
cor responding t o the orbital period of t he third b ody. 



IrwinI (119591 ) improved a method by IWoltjen (119221 ) for analyzing the long-term vari 



ation of the times of minima caused by a third body orbiting an eclipsing pair. Very 
use ful cora r nents and limitations were discussed by iFrieboes-Conde fc Herczegl (1l973l ) and 
by iMayer (Il99nh . Nowadays there a re several dozens of EBs, where the LITE is cer- 
tainly presented or supposed (see e.g. iBorkovits fc Hegedueslll996l . lAlbayrak et al.lll999l . 
Wolf et allliooi ). 

From the numerical point of view the method is a classical inverse problem. We have 
M measurements of the times of minima of the system at certain constant JDi with the 
individual uncertainties am,i- Our task is to find five parameters describing the orbit of the 
third body in the system: the period of the third body ps, the LITE semiamplitude A, the 
eccentricity e, the time of the periastron passage Tq, and the longitude of periastron uji2 for 
the binary on its orbit around the common barycenter. We have to compute simultaneously 
also two (or three) parameters of the eclipsing binary itself, namely its linear (or quadratic) 
ephemeris JDq and period P (and q for the quadratic one). Altogether, one has 7 (or 8) 
parameters to derive from the model fit of the minimum-time measurements 

{(JA,o-m,i)}i=i,A/ (p3,v4,e,To,u;i2, JL'o,^,?)- (2.1) 

The least-squares method and the simplex algorithm (see e.g. iKallrath fc Linnelll [l987h 
are used. The basic mathematic equations are the following. Compute the mean anomaly 
from the time of the measurement (the subscript i was omitted for the sake of brevity) 

M = 27r • (2.2) 

Then solve the Kepler equation M ^ E and convert the eccentric anomaly to the true 
anomaly of the third body in its orbit 



2 ■ arctan 



'1 + e E 

■ tan — 

1 - e 2 



After then, one can use the formula 
Ar ^ 



\J\ — cos^ 10^2 



1 



1 + e cos V 



sin(z/ + uj\2) + esina;i2 



(2.3) 



(2.4) 



to compute the magnitude of the LITE. Now it is possible to calculate the difference 
between the observed and calculated time of minimum 

{O - C) = JD - JDo - P ■ E - q ■ - At, (2.5) 

where E is the epoch of the JD according to the ephemeris JDq and P, and q is the 
quadratic term quotient. The resultant sum of normalized square residuals is 

Xhr.-t('-^t (2.6) 

Our task is to minimize this value and find the set of parameters (pa, A, e, Tq, cuu, JDq, P, q) 
describing the orbit. 

The weighting is provided by the uncertainties am,i- These values are obtained from 
the observations, or estimated as some typical uncertainty level for the certain kind of 
measurement provided by specific instrument. Another way is the following method. If we 
have information about the type of the observation (the method by which the measurement 
was obtained), we could use some weighting scheme Wi instead of uncertainties (e.g. Wi = 1 
for visual and Wi = 10 for photoelectric/CCD measurements) and solve the corresponding 
problem. Prom this solution we could find the uncertainties am,i simply as a differences 
between the observed and the predicted values. 

The third body in the system also causes variations in gamma velocities. If we know 
from different RV investigations and in different epochs of the system, we could see a 
variation of with a period corresponding to the orbital period of the third body. The 
variation could be described by 

Vj — K[cos{u + uu) + e cos 1012] , (2.7) 

where ir[km-s~^] is the amphtude of such variation and could be calculated from the LITE 
parameters A[d], psfyr], e and cuu from the equation 

X = (2.8) 

P3 V(l-e')(l-e2 cos2a;i2) 

But the basic limitation is very often the long period ps, which is usually too long to have 
reliable RV data for this analysis. 

LITE hypothesis could be also applied to other components in the system. The third 
body may cause LITE3 and another (fourth) component cause LITE4. The resultant total 
effect is then simply the sum of the two effects LITE = LITE3 + LITE4. One necessary 
condition has to be satisfied. The fourth component has to be more distant then the third 
one (and the third one distant from the EB pair). This is the main physical condition, 
because the method itself was derived with this condition by the use of von Zeipel's method 
to the three-body problem. 

If one wishes to include other effects, which could play a role in the problem, additional 
terms could be easily added. This means if one wants to describe a system with the LITE 
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caused by the third and the fourth component, also an apsidal motion is presented and 
mass transfer together, the variation in O — C diagram could be then described as 



(0-C) = JD-JDo-P-E-q-E^-{0- C)ute,3 - {O - C)utea - {O - C) 



apsid- 



One has to distinguish between primary and secondar y minima in apsidal motion term , 



which could be computed according to equations from lGimenez &: Garcia- Pelayd (119831 ) 



2.4 Astrometry 

Another method to study binaries and the properties of their orbits is astrometry. The 
number of visual binaries with astrometric orbits has grown, but complete phase coverage is 
often unavailable, due to the long orbital period involved. Thanks to precise interferometry 
the observable semimajor axes of astrometric binaries are still decreasing down to milli- 
and micro- arcseconds. On the other hand, one has to regret that no recent astrometric 
measurements of a wide pair of about 1" have been obtained for the systems mentioned 
below. Most of the astrometric observations were adopted from The Washing ton Double 



Star Catalogue WDS, see |http : / / ad . usno . navy. mil / wds / 1 (jMason et al.ll200ll ). The first 
astrometric observations of visual doubles are a few centuries old. Altogether there are 
about 2000 systems with their visual orbits known. 

The astrometric measurements of binaries consist of a series of measurements of po- 
sition angle 6i and separation pi secured at different times ti [i = 1,N). Sometimes also 
the errors of the individual data points are available. The weighting scheme is provided 
by using the uncertainties ae and dp of the individual observations. 

From the astrometric data one is trying to find the parameters of the relative orbit, 
defined by 7 parameters: period p^, angular semimajor axis a, inclination i, eccentricity e, 
longitude of the periastron tus for the third body on its orbit, the longitude of the ascending 
node Q, and the time of the periastron passage Tq. One has to solve the inverse problem 

{{ti, 6i, Pi, ae,i, (Jp^i)}i=i^N (a,p3,z,e,u;3,f2,To). (2.9) 

The least-squares method and the simplex algorithm were used. The basic mathematic 
equations are the following. Compute the mean anomaly from the time of the measure- 
ment, according to Eq. 12.21 solve the Kepler equation M ^ E and after then convert the 
eccentric anomaly to the true anomaly, according to Eq. 12.31 Compute the radius vector 
in arcseconds (the subscript i was omitted for the clarity) 

1 -e^ 

r = a-— , (2.10) 

1 + e cos V 

and from this equation one can compute the position on the sky, 9 and p, respectively: 



tan(^ — fi) = tan(z/ + cj) • cosi 

p = r ■ cos^u + u) sec{9 — Q) . 



(2.11) 
(2.12) 




Figure 2.1: A simplified description of the relative binary orbit on a plane of the sky. The true orbit is 
inclined against the plane of the sky (angle i). B denotes the barycenter of the whole system, 12 denotes 
the eclipsing binary pair and 3 the third component. The picture shows the moment, when the bodies 
are in the apocenters on their respective orbits. Also the semimajor axes of the third body and EB on 
the long orbit (ai2 and 03) are shown. 



Comparing this theoretical position on the sky with the observed ones 6*0 and po, one can 
calculate the sum of normalized residuals squared 



Xlstr 



N 
i=l 



9i — 6'( 



0,i 



Pi - PO,i 



a 



(2.13) 



foUowinglTorrea (j2004j ). With this Xastr using the simplex algorithm (see e.g. iKallrath fc Linnell 
19871 ) one can obtain a set of parameters (a, pa, 2, e, tua, fi, Tq) describing the astrometric 
orbit. 

At this place it is necessary to remark one useful comment. One has to distinguish 
between the two angles and oou- The parameter used in LITE analysis is uju, but in 
astrometry the quantity U3 = + tt is employed. In this thesis the angle u stands for 
the longitude of the periastron for the eclipsing binary, i.e. u = UJ12, and the subscripts 
will be omitted for clarity. 
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2.5 Combining the methods 

Our task is to combine the astrometry and the analysis of times of minima into one joint 
solution. If one has N astrometric and M minimum-time measurements, it is possible to 
merge them together and obtain a common set of parameters 



{U, di, Pi, (Tg^i, ap^i, JDi, a. 



m,i ) 



(A,p3,i,e,u;,fi,To, JDo,P,q). 



(2.14) 



This set of 10 parameters fully describes the orbit of the eclipsing binary around the 
common center of mass with the third unresolved component together with the ephemeris 
of the binary itself. 

One is also able to determine the mass of the third body and the semimajor axis of 
the wide system because the inclination is known and one can calculate the mass function 
of the wide orbit 



/(Ms 



[ai2 Sim) 



(Mgsini) 



(Ml + Ms + M3)2 pI 



173.15 -A 

Vl — cos^ cij 



[2.15) 



where ai2 stands for the semimajor axis of the binary orbit around the common center of 
mass and Mi,M2,M3 are the masses of the primar y, secondary, and tertiary component, 
respectively. For more details see e.g. Mayer 



The only difficulty which remains unclear is the connection between the angular semi- 
major axis a and the LITE amplitude A. The quantity ai2 could be derived from Eq. 12.151 
and with the masses of the individual components one is able to calculate also the value 
as, i.e. the semimajor axis of the third component around the barycenter of the system 



^3 — 0'12 ■ 



M 1 + M2 

M, ■ 



(2.16) 



The total mutual sep aration of the cornponen ts is atotai = ^12 + ^3 (see Fig J2.1l) . Using 
Hipparcos parallax tt (jPerryman fc ESAl 119971 ) one can obtain the distance D to the sys- 
tem. Now it is possible to enumerate the angular semimajor axis a as a function of D and 

O-total 

a = arcsin ^ _ (2.17) 

The w ay in which t he two different approaches were combined follows a similar ap- 
proach by iTorred (l2004l ). From the mathematical point of view both methods are analo- 
gous and there is an overlap of the parameters in both methods. Our task is to minimize 
the combined 

Xlomb — Xlstr + XLITEy (2.18) 

where xLtr ^i-nd Xlite sums of squares according to Eqs. 12.131 and 12.61 

There are some circumstances in which the values determined by the two methods 
are inconsistent. Especially when there are many more data points in one method than the 
other, the resultant would be much larger and as a consequence this method outweighs 
the other one. This problem could be eliminated using new uncertainties a instead of the 
old ones 



a 



M ■ Or 
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2.5.1 Error estimation 



The errors of the output parameters were calculated a ccord ing to the method described 
in Numerical Recipes, pages 684 - 694, see lPress et al.l (119861 ). Using a confidence limit of 
95%, one could calculate 

l-0.95 = r(z//2,AxV2), 

where F is the incomplete gamma function, u is the number of parameters fitted and Ax^ 
defines the boundary around the final solution. The area within this boundary is scanned 
and the maximum value of difference between actual parameter and final parameter 5a = 
{oi — ao) is taken as an error of the particular parameter. 



2.6 Distance determination 

Another task is determining the distance to these systems, which could be also done by 
combining the astrometry and the LITE analysis. As was mentioned above, from the 
LITE analysis it is possible to derive the quantities ai2sin2 and also Mssini, and from 
the astrometric analysis it is possible to compute the angle i and determine the semimajor 
axis and the third mass in absolute units. Thanks to these values the total semimajor 
axis atotai (see above) could be determined. Comparing the value atotai with the angular 
semimajor axis of the binary on the sky and leaving the parallax of the system as another 
free parameter, one could compute this value and derive the distance to the system. 

But this method is useful only in very special cases, where at least one period is 
covered with sufficient data points for both methods and both methods give us precise 
and comparable results. If the methods produce different results, the method is not useful 
for the distance determination. 



2.7 Limitations of the methods 

The methods themselves are very powerful and useful, but one has to take into consider- 
ation also some physical and observational limitations. 

If one considers only the LITE, the main observational limitation are the amplitude A 
and period p^. There are a few often-used methods to det ermine the time of minimurn , 



the most common being the Kwee-van Woerden method (IKwee fc van Woerderu Il956l ). 
Its main advantage is that one can compare the results and the individual errors of the 
measurements. It is suitable only for symmetric minima. 

For different types of eclipsing binaries we could reach different levels of accuracy of the 
times of minima. Using precise photoelectric and CCD detectors one is able to compute 
the time of minimum light with precision of about 1-10 seconds or less (~ 0.0001-0.00001 
day). This is only the theoretical value, attainable only if there are no clouds or moon, 
and if the observational conditions are very good. 

The practice is sometimes quite different. If one compares two observations of the 
same star and the same minimum time (after transformation to the heliocentric Julian 
time), one finds out that there could be a difference between them of order 0.005 day. 



I I 




Figure 2.2: In the left figure are the physical limits of amplitude of LITE as a function of period and 
mass. The different lines represent different masses of the components (for Mi = M2 = M^ ) . In th e right 
one the same for amplitude of systemic velocity variations. The figure adopted from iMayed ( 199Cl[ ). 



So there is a question about the accuracy of the method and the true error. Especially 
amateur astronomers have sometimes very precise measurements, but their results (the 
times of minima) are not very satisfactory. This could be due to shift of time on their 
computer (the delay could be from a few to tens of seconds!), or the wrong method used 
for determining the time of minimum. Another possible explanation is, that the time of 
exposure could be the beginning of the exposure (instead of the middle of exposure). 

Especially due to this reason there is a principal limit of the amplitude of LITE which 
could be reached. Amplitudes of LITE below 0.01 day are problematic, but in some systems 
even lower amplitudes are detectable (e.g. RT And or RZ Com). It is also necessary to 
take into consideration that the old measurements are not photoelectric, but visual or 
photographic, where the errors of the individual data points are much larger (and these 
errors are mostly not available for the analysis). 

Considering only LITE, plots of the relative limitations are shown in Fig. 12. 2[ One 
can see the dependence of the amplitude and the period on the masses of the individual 
components (only the case with equal masses is shown). 

If one uses the combined approach, it is necessary to discuss also the astrometric orbit 
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Figure 2.3: The principal limits of amplitude of astrometric orbit and amplitude of LITE in combined 
solution. The left figure was plotted with fixed LITE amplitude A = 0.03 d and the period = 20 yr. 
The right one with the same period and the amplitude of astrometric orbit a = 1 arcsec. The value of 
parallax was fixed, tt = 35 mas. The different curves represent the different masses of the EB, where the 
introduced masses are M12 = Mi + M2. 

and the accuracy of the astrometric measurements. A speckle interferometric techniques 
are very precise and could reach a few miliarcseconds (mas), but the older data are visual 
or photographic. In most cases the errors of individual data points are unavailable. 

Figure 12^31 shows the combined solution and the amplitudes in both methods. They are 
strongly dependent on the inclination i of the orbits. In the first one the LITE amplitude 
was fixed, while in the second one the astrometric-orbit amplitude was fixed. The basic 
properties of the plots in the diagrams are due to the geometry of the system. The edge-on 
orbit (i — > 90°) leads to the binary eclipses, while the face-on orbit {i — >• 0°) leads to the 
most pronounced astrometric variation. The parallax (the distance) of the system was 
also fixed. One has to consider also the limitations of the methods described above and 
the areas in the diagrams where these limits are satisfied. 
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Figure 2.4: The sum of squares of residuals as a function of number of iterations, for the case VW Cep 
Solution I. (see below). 



2.8 Numerics and the strategy to solve the problem 

The efficiency of tfie combined metfiod and tlie computing time required by tfie algoritlim 
strongly depend on tlie initial set of parameters and the input data. If nothing is known 
about the solution, one has to scan a wide range of parameters (eccentricity e from to 
1, and the angular parameters from 0° to 360°, etc.). 

The efficiency of the algorithm could be improved if the simplex is used repeatedly. 
It can happen that the simplex converges into a local minimum while the global one is 
far away. It is therefore advisable to run the algorithm again, with as large initial steps 
as in the previous run, but keeping the values of the parameters corresponding to the 
previously found minimum as one vertex. Repeating this strategy several times over the 
whole parameter space, one can judge whether the global minimum was found by checking 
whether the sum of squares of the residuals is still changing or not (see Fig. 12.41 and 12.51) . 

If one could guess the approximate values of the parameters, it is also recommended 
to use them. If one does not have any information about their values, the algorithm 
itself is able to find the appropriate ones, but these could be only numerical ones without 
any physical meaning. It is recommended to set the initial values of parameters as close 
as possible to the right values, because the algorithm itself will converge faster and the 
probability of the code producing incorrect parameters is lower. 

Another problem is the number of parameters used. The number of parameters strongly 
affects the computing time required. Regarding the classical LITE problem, one has 7 
parameters (2 from ephemeris and 5 from LITE), but using the quadratic term also one 
further parameter is necessary to evaluate. One can estimate the ephemeris parameters 
{JDq, P, q) in the first step and after then with fixed ephemeris calculate the parameters of 
LITE itself. But it is strongly recommended to compute all parameters together, because 
the LITE could also affect the ephemeris and change the values slightly. On the other 
hand sometimes it is better to fix the values of some parameters for the code to run faster. 
However, it is necessary to release all of the parameters for the final fit. 

If also another effect appears in the analysis, for example apsidal motion, it is necessary 
to estimate 3 additional parameters (cjq, e', cj). For the combined analysis, there could be 



J- I 




eccentricity 1 50 

omega(deg) 

Figure 2.5: The sum of squares of residuals as a function of cu and e for the third-body orbit. 

even more parameters. In the most complicated case (VW Cep, see section 1^^ . there 
were 14 parameters which had to be determined. The strategy which helps the code to 
converge faster to the minimum could be shown in the case HT Vir (see below section 
14. 5p . The astrometric orbit is known with high precision, so the astrometric parameters 
(a,p3, i, e, tu, f2, To) could be set and fixed for the first step. After the code returns the 
values of the ephemeris of the EB, one could release also the astrometric parameters and 
run the code once again with fitting all the parameters. This strategy saves a lot of 
unnecessary computing time. 

To conclude, using the strategy presented here and the combined method described 
above, one gets a satisfactory result after a large number of iterations. This number and 
computing time is strongly dependent on the input data set (the number of observed data) 
and the number of parameters fitted. For the case VW Cep with the largest data set (more 
than 1600 data points, see below) and also with the most parameters to fit (14 in total) 
one reaches the solution, when the sum of squares is not changing significantly, after circa 
100 000 simplex steps. This takes about one day on a computer with a 2 GHz processor. 
But it is only the illustrative example, the basic condition is the separation between the 
initial and the final parameters. 

There were done also a few tests of the code and its ability to find the same final 
parameters from different starting values of parameters. The results were satisfactory, 
and the code was able to find almost the same values, but it strongly depends on the 
separation between the initial and the final parameters. Also these tests confirmed the 
fact that setting the initial parameters as close to the right solution spares computing time 
and also makes the final solution more reliable. 
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2.9 The program 

The code itself was written in the Matlab language. It is available for download via the 
web pages http://sirrah.troja.mff.cuni.cz/~zasche/ The code is zipped together with the 
necessary routine files * . m, the sample input data files * . epo and * . dat and with the 
initial parameters file * . in. All of these files have to be copied into the same folder as the 
code itself. The sample (data files and also the code file) is for the HT Vir system. The 
code is easily modifiable for the user. A brief manual for the user is also available. 

The program was initially designed only for the use of this thesis, so the first version 
was not very user-friendly, but some modifications were done for the easier use. Some 
comment lines were also included in the file. 

In the first part of the code (lines 1 - 66) is the input. The files * . in, *.dat and 
* . epo include the input parameters, astrometry data and the minimum-time observations, 
respectively. The recommended data format is shown in the enclosed files. In the next 
lines (67 - 121) are some transformation rules and assigning the values to proper variables. 
Important here are lines 74-83, where is the input of masses Mi and M2 and also the 
parallax and its error. 

In the lines 122 - 172 is the relative astrometric orbit of the binary is plotted (before 
the computation) and saved as HTVir-before.eps. The next lines (173 - 299) are for 
plotting of the O — C diagram before the computation, saved as HTVirOC-before.eps. 

Lines 300 to 681 comprise the body of the code, the simplex algorithm. Its initialization 
is in line 327, where the user could choose which mode he wants to run. There are three 
possibilities: for only LITE, 1 for LITE together with the quadratic term (the mass 
transfer) and finally 2 for only the quadratic term (this possibility is not designed for the 
use of the combined approach). After then there are lines with the input parameters of 
the simplex algorithm (1.347 - 387). Here one could change the range of values of the 
individual parameters where the simplex works in the first run. The rest of the lines (to 
1.681) is the simplex itself and has not to be modified anyway. 

The lines from 683 to 1210 are for plotting the resulting — C diagrams and lines from 
1211 to 1268 for plotting the resulting orbit of the binary. In the lines from 1269 to 1322 
is the computation of the errors of the individual derived parameters. The lines from 1323 
to 1476 create the output files * . in and * .txt, where the output parameters are written. 

Generally, the code is ready to be used. The only modification could be the change of 
the names of the input and output files in the code (only find and replace the appropriate 
file names), and also the input masses and parallax. The main output file is *.txt, where 
are written all of the parameters with their respective errors. Also the derived quantities 
with the computed errors are included in the file. 

The computation process of the code is the following. Run the code * .m in the Matlab 
program. After two figures (0 — C diagram and astrometric orbit) appear on the screen, 
in the main Matlab command window the program asks the user whether he wants to 
compute also the quadratic term, or only the LITE. After confirmation which mode one 
wants to use, there appears another question about the number of iterations one wants 
to compute. It is up to user, but it is recommended to type only a few, because each of 
them could take some time. On the other hand more iterations means better precision of 
the result. After input these values on the screen will appear the resulting sum of square 



residuals of the problem, which is decreasing after each iteration. The sum of squares is 
divided into the two separate values, the first one from LITE and the second one from 
the astrometry. The most usual case is, that one of them is decreasing, while the other 
is increasing, but the sum of these two values has to decrease all the time. When the 
decreasing stops, the program terminates. A few plots will appear on the screen (these 
are saved as the * . eps files) and also the output file * . txt is saved to the same folder. 



Chapter 3 
Systems with LITE 



There were a lot of studies, where many echpsing binaries showing period variations of 
their period were analyzed. It was decided to present here only a few systems where LITE 
was not recognized until now, as well as systems where LITE was supposed to be present, 
but another (and better) solution for the variations in the O — C diagram has been found. 
This is not the crucial part of the thesis, and these systems are described only very briefly. 

Altogether there were about 130 systems analyzed during the 3-yrs PhD study for 
their period changes. In about one half of these systems the LITE was proposed as a 
hypothetical explanation. Regrettably, the analysis was done only on the basis of their 
times-of-minima observations and in most of the systems no other detailed analysis was 
performed. This is also evident in the set of binaries presented here in this chapter. 

A few of the analyzed systems we re selected f o r pub lication in various papers. Namely 



AD And, WY Pe r and V482 Per in IWolf et al.l mm . AR Aur, R CMa, FZ CMa and 



TX Her in IZaschd (120051) . 00 Aql, V338 Her, T LMi, RV Lyr, TW Lac and V396 Mon in 



Zasche et al.l (120061 ) . EW Lyr and IV Gas in lZaschd (120061 ) and XX Leo in lZasche fc Svoboda 
(j2006h . 



3.1 Individual systems under LITE analysis 



All of the selected systems are Algol-type EEs, and also semidetached one s. According 
to the recent paper on period changes in Algols by iHoffman et al.l (120061 ). there could 
be a connection between the spectral type of the secondary component and the nature 
of the period changes. Systems with spectral types of secondaries earlier than F5 show 
O — C variations, which could be cau sed by the r nagnetic activity c ycles and conyectiv e 
envelopes. This effect was discussed by lHalll (119891 ) . lApplegatd (119921 ). iLanza et al.l (119981 ). 
etc. The role of magnetic cycles on the period changes is discussed below, but due to lack 
of information about the systems such analysis is a difficult task. For some of the systems 
selected in this thesis the spectral types of secondaries are only known with a low confidence 
level, light-curve analysis is missing and spectroscopy has never been done. 

The L I TE an alysis of the systems presented below in this chapter was also published 
in IZaschel ( l2007h . 
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3.2 RY Aqr 



The eclipsing binary system RY Aqr (AN 125.1908, BD-11 5 574, HP 2030 69) is an Algol- 
type EB. It was classified as A3 spectral type (Simbad), A8 flPoDPeij 19891 ). or most likely 
as a late A /early F main sequence star and an early K subgiant fiHeltlll987l ). It is a double- 
lined SB (jPopperlll989l ) with an orbital period of about 2 days. Its relative brightness is 
about 8.9 mag in V filter. 



The variability of RY Aqr w as disc o vered by iLeavitt fc Pickering (119081 ) and the 
ephemeris was firstly derived by IZinned (119131 ). Since then a lot of times of minima 
observations were obtained. 



The most detailed analysis was performed by iHeltl 119871 on the basis of the uvh ii pho- 
toelectric photometry together with the radial velocity measurements by iPoppen This 
analysis (using both WINK and Wilson-Devinney programs) results in a set of parame- 
ters, which reveals the nature of the system. The hotter primary has unusually low mass 
(about 1.27 Mq), while the secondary K subgiant (0.26 Mq) has undergone a mass-loss 
from its initial mass of about 1.9 M0. There was observed also an intrinsic photometric 
variability which could be caused by the surface activity of the secondary (the period of 
this variability is close to the orbital period of the system). RY Aqr is also a member of 
a visual binary HU 86, but the distant body probably doe s not relate to t he systern. 

There were performed a few period studies of RY Aqr (jBaldwinlll973 . lMallamalll980l ) 
and the most recent one by iHeltl ( 119871 ) , who suggested that also the LITE could play a 
role in this system. Her supposed period of such variation (about 70 years) is different 
from the present one. This new analysis is based on a set of 178 times of minima published 
in the literature. Resultant O — C diagram is plotted in Fig. 13.11 and the parameters of 
the third-body orbit are in Table 13. 1[ 

The minimal mass of the third body (Ms^mm = ^3 ■ sin is if = 90°) was calculated 
according to the mass function derived from the LITE hypothesis and total mass of the 
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Figure 3.1: An O — C diagram of RY Aqr. Tlie individual observations are shown as dots (primary) 
and open circles (secondary), the small ones for visual and the large ones for CCD and photoelectric 
observations, bigger the symbol, bigger the weight. The curve represents the predicted LITE variation. 
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Table 3.1: The final results: RY Aqr. 



Parameter 


JDo 


P 


P3 


Tq 


UJ 


e 


A 




Al'i rri-i'n 

'J ^ 1 1 Lot L 


Unit 


[HJD] 


[day] 


[yr] 


[HJD] 


[deg] 




[day] 


[Ms] 


[Mo] 


Value 


2440824.351 


1.9665990 


105 


2442300 


88 


0.35 


0.070 


0.165 


1.02 


Error 


±0.004 


±0.0000013 


±4 


±1600 


±7 


±0.05 


±0.002 


±0.006 


±0.03 



eclipsing pair M12 = (1.27 + 0.26) Mq, see Table 13.11 Computed minimal mass about 
1 Mq is rather high and such a star on the main sequence will be evident bec ause of its 
lumin osity. The mass M3 „^j„ leads to the spectral type around G3 (according to lHarmanec 
19881 ). and such a star would be more luminous tha n the secondary component of the EB 



pair. The third light in the light-curve solution by (119871 ) was only estimated. The 
value about a few percent (from 2.6% in u to 5.9% in y) was adopted, but not derived. 
This value indicates roughly the same third mass, as was cal culated from o ur analysis. 

The radial velocities were analyzed precisely only once ( Popper 19891 ). therefore no 
changes in systemic velocities are available. The systemic velocity of RY Aqr was derived 
to be about —60 km ■ s~^, which could be caused by the motion around the common center 
of mass with the third component. Precise spectroscopy would probably detect the third 
body in the spectrum of the system. 

The star has not been m easured by Hipparcos, but the distance was derived from the 
photometry, see iHeltl (119871 ). The value (180 ± 10) pc leads to the predicted angular sepa- 
ration of the third component a = (169 ± 10) mas, which was calculated using assumption 
^3 = 90°. The predicted magnitude difference is about 3 mag. The companion with such 
a distance and magnitude difference is detectable with the modern stellar interferometers. 
Only further times of minima, precise spectroscopic and photometric analysis will reveal 
the nature of the system. 



3.3 BF CMi 



BF CMi is one of the neglected eclipsing binaries, which have been observed only a few 
times and only very limited knowledge about it is available. Its p eriod is about 1.18 days 
and th e relative brightness about 10.3 mag in V filter. According to lSvechnikov fc Kuznetsova 
(119901 ) the star has spectral type A5+K0IV (based only on photometric indices), mass ratio 
0.3, orbital inclination 79° and is a semidetached one . 

Its variab ility was discove r ed by Huruhatal ( 1979 ) and its designation as BF CMi was 



presented by iKholopov et al.l (Il98ll ). But the period is still questionable. There were two 
unsuccessf ul attemp t s to o bserve the secondary minima - on 2 March 2006 and 4 April 
2007. Also iBertholdl (jl98ll ) noted that no secondary minimum is observable. It means the 



Table 3.2: The final results: BF CMi. 



Parameter 


JDo 


P 


P3 


To 


a; 


e 


A 




-^^3, mm 


Unit 


[HJD] 


[day] 


[yr] 


[HJD] 


[deg] 




[day] 


[Mq] 


[Mq] 


Value 


2450789.610 


1.1806791 


46.3 


2447300 


170 


0.79 


0.040 


0.39 


2.1 


Error 


±0.016 


±0.0000026 


±1.2 


±200 


±21 


±0.10 


±0.018 


±0.16 


±0.9 
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Figure 3.2: An O — C diagram of BF CMi. The description is the same as in the previous figure. 



period could be 2 times longer, about 2.36 days. Two new primary minima were observed 
and kindly sent by L. Smelcer from Valasske Meziffcf observatory. 

There were 39 times of minima collected from the published literature, resulting in 
an O — C diagram shown in Fig J3.2[ The parameters of LITE are in Table 13 .21 As is 
evident from the abrupt period jump near 1990, the eccentricity of the orbit should be 
rather high. Note also quite large scatter of recent photoelectric and CCD observations 
since 1990, which is much larger, than one could expect from these kind of measurements. 

If one assume the masses of the individual components Mi + = (1.9 + O.9)M0, the 
minimal mass of the predicted third component is about 2.IM0, which is rather high value 
and would dominate in the system. This hypothesis could not be proved until the detailed 
analysis of the system is performed. Regrettably, neither photometry nor spectroscopy 
was carried out. The abrupt changes in period could be also caused by the two period 
jumps near 1988 and 1991, produced by some mass-transfer phenomena in the system. 



3.4 RW Cap 



The eclipsing binary RW Cap (AN 21.1910, BD-18 5641, HD 192900) is a system with an 
orbital period of about 3.4 days. Its spectral type was classified as A3+A4 (according to 
Budding||l984l ) and its apparent brightness is about 10.3 mag in V filter. The depth of 
prima ry minimum is about 1.2 mag. Stromgren photometry of the system by I Wolf fc Kern 
( 119831 ) agrees with its spectral type. 

Its photometric variability was discovered by Pickering ( 191ol ) and until now there were 
52 observations of times of minima obtained. IZessevichl (119571 ) c ollected all avai lable min- 
ima and proposed an abrupt period jump near 1920. After then iKreinerl (119711 ) compiled 
large set of times of minima and already in this paper is evident that there could be some 
periodic variation in O — C diagram. Regrettably during the last two decades only a few 
times of minima were obtained, so the LITE hypothesis is still not very conclusive. Two 
of last three data points in the O — C diagram are only a mean values from the automated 
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Table 3.3: The final results: RW Cap. 



Parameter 


JDa 


P 


P3 


To 


UJ 


e 


A 


<? 


f{M3) 




Unit 


[HJD] 


[day] 


[yr] 


[HJD] 


[deg] 




[day] 


[10-10 day] 


[Mo] 


[Mq] 


Value 


2435750.857 


3.3923745 


80.1 


2440900 





0.23 


0.130 


87.4 


1.91 


5.9 


Error 


±0.014 


±0.0000055 


±4.6 


±2800 


±42 


±0.16 


±0.011 


0.2 


±0.26 


±1.0 



surveys ROTSE (see lAkerlof et allbnnnl ) and ASAS (see lPojmanskilll997l ). But the last 



one was obtained and kindly send by A.Liakos from Athens university in July 2007. 

For the times of minima analysis the LITE and the quadratic term was used. This 
means, during the computation process, altogether 8 parameters {JDq, P, q,p3, A, To,a;, e) 
were adjusted, resulting in a set of parameters written in Table I3.3[ Because the system 
is semidetached, the mass-transfer hypothesis could play a role. The quadratic term 
coefficient q = (87.4 ± 0.2) ■ 10^^° day leads to the period change about 1.88 • 10^® day/yr. 
From this value the conservative mass transfer rate could be derived M = 9.4-10"'^ Mq/jt. 
Regrettably, these values cannot be proved by some other independent method. 
Applying the LITE hypothesis to the system one gets the 80-years variation 
but as is evident 



l^see 

Fig j3.3p . but as is evident from Table 13.31 the orbit is still not very well-defined and 
the errors of the individu al parameters are large. Res ulting mass function is quite high 
and using Mi2 = 4.59 My flBrancewicz fc DworaHll980h one gets the minimal mass of the 
third body of about 5.91VI0. Such a body would be dominant in the light-curve solution as 
well as in the spectroscopic analysis. Unfortunately there were no such analysis performed. 
Another explanation is that the third component is also a binary. 



3.5 TY Cap 

The next eclipsing binary with period changes is TY Cap (AN 243.1932, BD-13 5664, 
HD 194168). It is an Algol-type EB with apparent brightness about 10.3 in V filter and 
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Figure 3.3: An O — C diagram of RW Cap. The description is the same as in Fig l3.H the dash-dotted 
line represents the quadratic term. 
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Figure 3.4: An O - C diagram of TY Cap. 



spectral type classified as A2/3V (according to lHalbedellll984j ) . Its orbital period is about 
1.4 days. 

Its photometric variability as well as its Algol-type were discovered by iHoffmeister 
( 119331 ). Altogether 96 times of minimum light were carried out, only 5 data points were 
neglected due to their large scatter. The O — C plot is in Fig J3.4l the curve represents the 
least-square fit with the LITE parameters written in Table 13. 4[ 

Due to missing det ailed analysis of the sy s tem, o ur knowledge about TY Cap is only 
limited. According to iBrancewicz fc DworakI (119801 ) the total mass of the EB system is 
Ml + M2 = (2.5 + 2.O6)M0. With this mass and with the parameters of the LITE from 
Table [33] one can calculate the minimal mass of the third body, which results in 2.18 Mq. 
Unfortunately there is no spectroscopic, as well as no light-curve analysis and the star was 
not measured by Hipparcos, so the angular separation of the third component also cannot 
be derived. 

The eccentricity of the LITE solution is quite high, and could be even higher, but due 
to lack of data points near the periastron passage this could not be proved. The next 
periastron passage will occur about 2035. Generally the third-body orbit is not covered 
sufficiently and the parameters of this predicted body have to be derived by some other 
independent method. 



Table 3.4: The final results: TY Cap. 



Parameter 


JDo 


P 


Pa 


To 




e 


A 


/(M3) 


-^^3, mm 


Unit 


[HJD] 


[day] 


[yr] 


[HJD] 


[deg] 




[day] 


[Mq] 


[Mq] 


Value 


2444793.489 


1.4234574 


70.4 


2439600 


147 


0.79 


0.045 


0.23 


2.2 


Error 


±0.006 


±0.0000009 


±8.6 


±1200 


±29 


±0.12 


±0.005 


±0.02 


±0.2 
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Figure 3.5: An O — C diagram of SS Get. The description is the same as in the previous figures. 



3.6 SS Get 

SS Get (BD+01 491, HD 17513) is a semidetached Algol-type EB with an apparent bright- 
ness of about 9.4 ni ag in V filter and spectral type classified as A0+K3III (according to 
Budding et al.ll2004j ). Orbital perio d is about 3 days. 

Its variabi lity was discovered by iHoffmeisterl (1l934al ). The most detailed analysis was 
performed by lNarasaki fc Etzell (119941 ) on the basis of their BVRI photoelectric photom- 
etry and spectroscopy. This study results in a semidetached system, with no peculiarities 
in the light curve. Mi = 2.15M0, M2 = O.6M0, and the spectral observations indicate that 
no circums tellar matter i s pres ented in the system. On the other hand the spectroscopic 
analysis by IVesper et al.l (120011 ) presents evidence for mass transfer in the binary because 
of the behavior of the Ha emission. The mass transfer was not recognized in our current 
analysis of the times of minima. 

Since its discovery there were 111 minima observations obtained, but only 95 were 
used for the period analysis, because of the large scatter of the first ones from the 1930's. 
There could be a period jump near 1950, this is the reason why on ly the recent data were 
analyzed. The period changes were firstly noted by iKreinerl (jl97ll ). but only with a small 
set of times of minima, displaying the steady increase. The final LITE curve is in Fig. 13.51 
and the resultant parameters in Table 13.51 

Using the mass of SS Get derived by lNarasaki fc Etzell (119941 ) Mu = 2.75M0, one gets 



the minimal mass of the third component about O.72M0. The predicted value of the third 
light is only about 1 % and the third component would be also similar to the secondary 



Table 3.5: The final results: SS Get. 



Parameter 


JDo 


P 


P3 


To 




e 


A 


/(M3) 


-^^3, mm 


Unit 


[HJD] 


[day] 


[yr] 


[HJD] 


[deg] 




[day] 


[Mo] 


[Mo] 


Value 


2442451.330 


2.9739737 


20.4 


2449500 


304 


0.21 


0.014 


0.031 


0.75 


Error 


±0.002 


±0.0000007 


±0.5 


±900 


±73 


±0.09 


±0.002 


±0.002 


±0.13 
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Figure 3.6: An O — C diagram of TY Del. The description is the same as in the previous figures. 

component in the spectra. Detection of the third body in the hght curv e as well as in the 
spectr a is a difficult task. The distance to the system was derived by iNarasaki fc Etzel 
(jl994j ). resulting in c? = 486 pc. According to this value and derived parameters of the third 
body, one could calculate the predicted angular separation of the third component about 
only 23 mas and magnitude difference about 5 mag. Such a large magnitude difference 
and low separation is hardly detectable, and the third body remains undetectable also by 
interferometry. 



3.7 TY Del 



Another EB with apparently variable period is TY Del (AN 141.1935, BD+12 4539), spec- 
trum classified as B9+G0IV ( iHoffman et al.ll2006l ) and relative brightness about 10.1 mag 
in V ffiter. There is a consensus about the spe ctral types of the comp o nents of TY Del, 
but there is a difference between the masses. Brancewicz &: Dworak fll980h . and after 
then also Buddind ( 1984) and Budding et al. (2004) have presented the masses Mi = 5M0, 



Ma = 2M0, while ISvechnikov fc Kuznetsoval (119901 ) presented Mi = 2.8Mo, Ms = O.84M0, 



what is more likely to the proposed spectral types. 



The star was discovered to be a variable by iHoffmeisted (119351 ). Th ere was only on e 
attempt to observe the whole light curve of TY Del photoelectrically by iFaulkned (119831 ) , 
unfortunately only about half of the cu rve was obs erved. No analysis of these data was 
carried out. The star was also studied by lCookl (119931 ) on the basis of his visual observations 



Table 3.6: The final results: TY Del. 



Parameter 




P 




7i) 




e 


A 


/(M3) 




Unit 


[HJD] 


[day] 


[yr] 


[HJD] 


[deg] 




[day] 


[Mo] 


[Mo] 


Value 


2442959.471 


1.1911264 


64.9 


2449200 


38 


0.22 


0.027 


0.025 


0.79 


Error 


±0.001 


±0.0000002 


±2.3 


±1000 


±18 


±0.06 


±0.002 


±0.001 


±0.07 
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for the fong-time scale intrinsic variations, but the re sults are not very co nclusive. 



The spectroscopic observations in Ha were done by lVesper et al.l (l200ll ). They conclude 
that there is no activity in Ha and no evidence for the mass transfer structures was found 
in this system. 

Altogether 370 times of minima were collected, from which only 5 were omitted due 
to their large scatter. One period of the third body is already sufficiently covered by data 
points, see Fig. 13.61 but further observations are still needed. The last one data point was 
observed at Ondfejov observatory. From the LITE parameters (see Table and with the 
approximate m asses of the individual cornpone nts of the eclipsing binary Mi = 2.8 Mq, 
Ms = 0.84 M0 f jSvechnikov fc Kuznetso^^ll990h one is able to derive the minimal mass of 



the third component M^^min = 0.67 Mq. Due to lack of any other observations also this 
hypothesis cannot be proved. The spectral types and masses were derived only on the basis 
of the photometric indices and are not very conclusive. The spectroscopic analysis, as well 
as the analysis of the light curve of the system is needed, but the third light is undetectable 
in the light-curve solution. Regrettably the star was not measured by Hipparcos, so the 
distance is not known and one cannot derive the predicted angular separation of the third 
component. As one can see, there is some additional variation besides LITE, which is not 
strictly periodic, see Section 13.141 for details. 



3.8 RR Dra 



Another eclipsing binary which exhibits apparent period changes is RR Dra (AN 188.1904). 
It is an Algol-type sem idetached binary, relative bright ness ab out 9.8 mag in V, spectrum 
classified as A2-I-G8IV (ISvechnikov fc Kuznetsovalll990l). whilelYoon et al.l (Il994r) classified 
a little bit later spectral type of secondary KO. ISvechnikov &: Kuznetsoval (1l990l ) presented 
the masses Mi = 2. ISM© and M2 = O.6M0. The primary minimum is very deep, about 
3.5 mag and the orbital period is about 2.8 days. 




Figure 3.7: An O — C diagram of RR Dra. For the plot where the quadratic term was subtracted see 
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Table 3.7: The final results: RR Dra. 





T n 


p 


P3 


Jo 


UJ 




A 

/\ 


9 


Jv^hj 




Unit 


[HJD] 


[day] 


[yr] 


[HJD] 


[deg] 




[day] 


[10-1" day] 


[Mq] 


[Mo] 


Value 


2434913.728 


2.8312140 


84.3 


2450100 


110 


0.50 


0.073 


-126.2 


0.300 


1.85 


Error 


±0.022 


±0.0000053 


±0.6 


±400 


±4 


±0.03 


±0.002 


0.2 


±0.002 


±0.09 



The star was discovered to be a variable by IVliss ICeraskil (119051 ). The minimum is so 
deep that also visual observers could provide reliable observations. That is the reason why 
most of the collected times of minima are the visual ones (193 out of 219). iKreinerl ( 1l97ll ) 
collected all available minima for the period analysis. The long-time increase of the period 
is evident from his O — C diagram (due to the mass transfer bet ween the co r npone nts?). 
The most recent period study of this system was performed by iQian et al.l (120021 ) . who 
considered (besides the mass transfer) the abrupt period jumps - altogether 8 jumps were 
introduced to describe the O — C diagram in detail. Almost the same goodness of fit could 
be reached by applying the LITE hypothesis besides the mass transfer. 

Altogether 219 times of minima were used for the analysis. One new observation of 
minimum was obtained by IvI.Wolf at Ondfejov observatory. The O — C plot is in Fig J3.7l 
where LITE and the quadratic term were plotted together. In the next figure only LITE 
is shown, see Fig J3.8[ As one can see, the period increase is very rapid, and the amplitude 
of LITE is still quite high. This leads to the relatively high mass function, which results in 
high predicted minimal mass M^^rnin = I.SSIvIq. This is larger than the secondary and in 
the light-curve solution as well as in the spectrum will be probably observable. Regrettably 
no such analysis was performed. 

The quadratic term coefficient q = (126.2 ± 0.1) ■ 10"^'^ day leads to the period change 
about 3.26 ■ 10~^ day/yr. From this value the conservative mass transfer rate could be 
derived M = 3.5 ■ 10"'^ Mq/jt. This relatively high value of mass transfer rate arises 
from the very rapid period change, which was attributed to the q uadratic ephemer i s. The 
spectroscopic observations during the primary eclipse made by iKaitchuck et al.l (119851 ) 
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Figure 3.8: An O — C diagram of RR Dra after subtraction of the quadratic term. 
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Figure 3.9: An O - C diagram of TZ Eri. 



indicate the possible presence of a transient accretion disc in the system. The presence of 
such a disc also supports the hypothesis of mass transfer in the system. 

For the estimation of the angular separation of the third body and the astrometric 
confirmation of the LITE hypothesis the distance to the system has to be known. The 
star was not measured by Hipparcos and the distance is not known precisely. The only 
information about the distance is from iKharchenkol (j2nnih . where is introduced a surpris- 
ingly inaccurate value of the parallax vr = (0.40 it 11.50) mas. Distance with such a large 
error is useless for the estimation of the angular separation of the predicted component. 



3.9 TZ Eri 



The system TZ Eri (AN 40.1929, BD-06 880) is an EB with an orbital period of about 2.6 
days, apparent brightness of about 9.7 mag in V filter. It has a deep primary minimum 
(about 2.8 mag), so the visual obse rvations could be a lso reliable. 

Its variability was discovered by iHoffmeisterl (1l929l ). who also r ecognize d the s ystem 
to be an Algol-type. The spectral type was first classified by Miss ICannonl fll934h as r . 
Later the spec t rum was re-classified as A5/6 V (primary) and KO/1 III (secondary) by 



Barblan et al.l (Il998l ). In this later paper the light-curve observations in the Geneva 7- 
colour photometric system were analyzed together with the radial- velocity curves of both 
components. Wilson-Devinney code was used, resulting in a set of parameters describing 
both components. For our analysis are the most important the masses, M\ = 1.97 M0 and 
M2 = 0.37 M0. 



Table 3.8: The final results: TZ Eri. 



Parameter 


JDo 


P 


P3 


To 


UJ 


e 


A 


<? 


/(M3) 




Unit 


[HJD] 


[day] 


[yr] 


[HJD] 


[deg] 




[day] 


[10-1" day] 


[Mo] 


[Mq] 


Value 


2446109.730 


2.6061129 


48.8 


2451100 





0.01 


0.042 


-18.0 


0.165 


1.3 


Error 


±0.009 


±0.0000034 


±6.8 


±2300 


±40 


±0.10 


±0.014 


0.2 


±0.013 


±0.1 
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T here were also several studies about the presence of the accretion disc in the system 
(e.g. iKaitchuck fc Honevcuttlll982l . lKaitchuk fc Parklll988l . IVesper et aLllioOlh . This disc 
as well as mass transfer from the secondary to the primary is in agreement with our 
result about the increasing orbital period (see below). The system was a lso included in 
the sample of Algol- type binaries with radio emission ( jUmana et al.lll998h . The star was 
also investig ated according to the p ossible connection between the orbital and pulsational 
periods, see ISoydugan et al.l (120061 ). 

The analysis of the long-term period changes was done with a set of 108 observations 
(mostly the visual ones). The resultant O — C diagram is in Fig l3.9l and the parameters 
of the predicted LITE are in Table [331 The minimal mass of the third corn ponent results 
in Ms^mm = 1-3 Mq, or the spectral type F6 (according to lHarmanedll988l ). Such a body 
could be evident in the light-curve solution as well as in the spectra of TZ Er i. Regrettably, 



there was no attempt to detect such a body during the detailed analysis by iBarblan et al. 
( 119981 ). The long-term period increase is due to the mass transfer from the secondary, 
with the conservative mass-transfer rate M = 6.2 ■ lO~^M0/yr.^ 



Despite the fact the star was not observed by Hipparcos, iBarblan et al.l (119981 ) esti- 
mated the photometric distance to d = (270 ± 12) pc. Assuming the coplanar orbit of 
the third component, then M3 = M^^rnin and one could calculate the predicted angular 
separation of the third body to a = 77 mas and magnitude difference about 1.7 mag. Such 
a component is detectable with the modern stellar interferometers. 



3.10 RV Per 



The system RV Per (AN 61.1905, BD+33 805, HD 279552) is an EE with an orbital period 
of abo ut 2 days, spectral type classified as A2+G7IV, according to lSvechnikov &: Kuznetsova 
( I1990I ). Also this star shows deep primary minim um, about 2.4 mag. 

The star was discovered to be a variable by iBlazkol (119071 ). Since then only a few 
papers on this star were published, so our knowledge about the system is very hmited. 




Figure 3.10: An O — C diagram of RV Per. 
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Table 3.9: The final results: RV Per. 



Parameter 


JDo 


P 


P3 


To 




e 


A 


/(M3) 




Unit 


[HJD] 


[day] 


[yr] 


[HJD] 


[deg] 




[day] 


[Mq] 


[Mq] 


Value 


2442046.920 


1.9734888 


99.8 


2431600 


210 


0.79 


0.017 


0.007 


0.47 


Error 


±0.002 


±0.0000003 


±18.8 


±7200 


±44 


±0.46 


±0.003 


±0.002 


±0.07 



According to iBrancewicz fc DworakI (jl980l ) the masses of the individual components are 
Ml = 3.04 M0 and = 0.46 Mrr^. The attempts to prove the existence of the accretion 
disc in the system ( iKaitchuck et aLlll985l ) were not success ful, as well as th e presence of 
the pulsating component in the system was no t confir n aed (jKim et al.ll2003l ). 

The period changes were first studied by IWoodI (Il950l ). but the data set was not 
sufficient to do any satisfactory conclusions. In the present thesis the data set consists 
of 146 times-of- minima observations (see Fig iS.lOp . One new minimum was observed by 
M.Zejda. The parameters of LITE are in Table [3l9l As one can see, the period of the third 
body is not covered by observations yet and the errors of the individual parameters are 
high. The value of eccentricity could be even higher, but there are no observations near 
the periastron passage, and the next one is predicted to occur near 2040. Only further 
observations would confirm or reject the third-body hypothesis. 



3.11 UZ Sge 



The Algol-type EB system UZ Sge (AN 435 .1936) has an orbital period of a bout 2.2 days 
and spectral type classified as A3V+G0IV fISvechnikov fc Kuznetsov^ll990h 

Its photometric variability was discovered by lGuthnick fc Schnellerl ( 119391 ). Since then 
there was no attempt to do any detailed analysis, neither th e photona e tric ri or the spectro- 
scopic one. The only spectroscopic observation was done by iHalbedell (Il984j ) for derivation 
of the spectral type of primary component. 

Altogether 122 measurements of times of minima were found in literature, but 14 
observations were neglected. Four new observations of minima were obtained (two of 
them by L. Smelcer, one by M.Wolf and one by author ) . If t he masses of the individual 
components were taken from ISvechnikov fc Kuznetsoval (119901 ). Mi = 2.05 Mq and M2 = 
0.29 M0, then the minimal mass of the third component results in Ms ^m = 0.65 Mq. But 
due to absence of any detailed analysis of this system, this value cannot be proved. 



Table 3.10: The final results: UZ Sge. 
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Figure 3.11: An O - C diagram of UZ Sge. 



3.12 BO Vul 

The last EB system in this thesis which shows long-term period changes is BO Vul 
(AN 125.1935, HD 345287). It is an eclipsing binary with an orbital period of about 
1.9 days, apparent brig htness 10 mag, depth of primar y minimum about 1.6 mag and 
spectral type FO+GOIV f Svechnikov &: Kuznetsova 1990l). 

The star was o bserved to be a variable by iHoffmeister fll935h and the first ephemeris 
were presented by iGuthnick fc Pragerl (119361) . The first observation of the whole light 
curve was carried out by iNassaul (jl939l ). where also a brief analysis of the system was 
presented. Since then there was no detailed analysis of t he syste n i perf ormed. 

The c hange s of its period were firstly mentioned by lAhnertl (119731 ). After then also 



BaldwinI (119961 ) published new revised elements for BO Vul, but without any interpretation 



of the period changes. 

Altogether 390 times of minima were collected, but only 360 were used for this analysis. 
One new minimum was observed by M.Wolf. For the fin al results see Fig ]3.12l and Table 
13. Ill If the masses from ISvechnikov fc Kuznetsova! (1l990l ) were taken. Mi = 1.45 Mq and 
M2 = 0.64 Mq, then the minimal mass of the third component results in 0.73 Mq. But 
as in the previous cases, there is no detailed analysis, which could prove this result. The 
quadratic term leads to the conservative mass-transfer rate of about M = 1.3-10^^ Mq/jt. 

One can also see some additional non-periodic changes, which are evident since 1970's 
and which could not be described by applying only LITE and the quadratic term. The 
amplitude of these variations is smaller than the amplitude of LITE, but one cannot doubt 
about their presence nowadays. These could be caused by the abrupt period changes, or 



Table 3.11: The final results: BO Vul. 
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Figure 3.12: An O - C diagram of BO Vul. 



due to magnetic activity cycles presented in the system. See the next section for a brief 
analysis. 

According to iBrancewicz fc Dworak the distance to the system is about 233 pc, 

but the star was not measured by Hipparcos, so this value is not very reliable (there is 
no information about the error of this value). From the distance one could estimate the 
predicted angular separation of the third component, resulting in 63 mas, and magnitude 
difference about 3.4 mag. Such body is perhaps marginally detectable by the modern 
stellar interferometers. 



3.13 Alternative explanation 



One can also see additional non-periodic variations in some of the O — C diagrams, which 
could not be described by applying only the LITE hypothesis. In Figs. 13.131 there are 
shown four cases with the most evident variations. The amplitudes of these variations are 
usually about 10 minutes in the — C diagram and are not strictly periodic (the "periods" 
are from 5 to 20 years). This could be caused by the presence of stellar convection zones 
and magnetic activity cycles in an agreem en t with so-called Appl egate's mechanism, see 
e.g. lApplegatel (jl992l ). iLanza et all (jl998h . or lHoffman et all (j2006h . The effect could play 
a role, because the spectral type s of th e secondary components in most of the systems are 



later than F5 (see I Zavala et al.l (120021 ) for a detailed analysis) . This explanation would 



clarify the non-periodicity and the changes in amplitude of such variation, as well as why 
in some binaries this phenomena is presented, while in some others not. 

Due to missing information about the properties of the eclipsing components in most 
of the systems one also cannot estimate the predicted variation of the quadruple moment 
AQ, which causes the period variations. This value could be computed from the equation 



AP 



-9 



AQ 



where P is the orbital period of the system, M is a mass of the star and a is the separation 
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Figure 3.13: The O — C diagrams of four systems after subtraction of the LITE and the quadratic term 
(cases RR Dra and BO Vul). The additional variation is clearly visible. 



between the components, see iLanza fc Rodond (120021 ). AP is the amphtude of the period 
oscillation and could be computed from the LITE parameters from equation 



AP = A- 



2[l-cos( )], 



P3 
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The typical values of AQ are of the order of 10 
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Due to missing light-curve and radial-velocity curve analysis, the value a is missing. 
For the few cases where this value is known (RY Aqr, SS Cet and TZ Eri) only RY Aqr 
does not satisfy the condition about the limits for AQ (being about 10 times lower). This 
result does not indicate that the magnetic activity cycles are not presented in this system, 
but only the fact that this effect cannot be used as an alternative explanation of the period 
changes. The effect could be present in addition to the light-time effect and describe the 
non-periodic variations (shown in Figs. 13.131 after subtraction of LITE). 

To conclude, for the better description of the observed period variations of these sys- 
tems, the magnetic activity cycles could be presented together with the LITE. On the other 
hand one has to take into consideration that the spectral type s of most of these binaries 



were derived only on the basis of their photometric indices (jSvechnikov fc Kuznetsova 
1990) and are not very reliable. 



3.14 Brief summary 

Eleven Algol-type semidetached eclipsing binaries were analyzed for the presence of LITE 
on the basis of their O — C diagram analysis and the times-of-minima variations. A few 
new observations of these systems were obtained. All of the systems above show apparent 
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changes of their orbital periods, which could be explained as a result of orbiting the EB 
around the common center of mass with the third component. 

Such a variation usually has a period on the order of decades. The light-time effect was 
applied as a main cause of these changes (as one can see from Figs j3.lH3.12l) . In four cases 
(RW Cap, RR Dra, TZ Eri and BO Vul) also the quadratic term in the light elements was 
used. This could be described as a mass transfer between the components, which could 
play a role, because all of the systems are semidetached ones. In some cases also the proof 
of presence of mass-transfer structures or accretion discs were revealed by spectroscopy. 

Regrettably, in most of the systems no detailed analysis (neither the photometric nor 
the spectroscopic) was carried out. The spectral types and the masses of the individual 
components in most of the systems are only approximate, so the parameters of the pre- 
dicted third bodies are also affected by relatively large errors. Due to missing information 
about the distances to most of these binaries also the prediction about the angular sep- 
aration could not be done. As one can see, only further detailed photometric, as well 
as spectroscopic and interferometric analysis would reveal the nature of the system and 
confirm or refute the third-body hypothesis. 



Chapter 4 

Systems with combined LITE and 
astrometry 

The crucial part of this thesis was the analysis of the systems, where the EB pair is a 
component of spatially resolved binary. Despite increasing number of the visual as well as 
eclipsing binaries, the intersection of these two sets is still only very limited. 

Finding appropriate candidates for this analysis turned out to be quite difficult. One of 
the problems was the data set and its quality. Such systems have to satisfy the following 
adopted conditions: 1. More than 10 times of minima and more than 10 astrometric 
observations are available. 2. The observed range of the position angle in the astrometric 
measurements is larger than 10°. The limit for the number of data points was accepted 
because of the number of parameters, which have to be found. There were 5 parameters 
for LITE and 2 ephemeris, or 7 for astrometry, it means at least 7 data points (in both 
methods) is the absolute minimum for the analysis. And the limit for the range of data 
points was the accuracy of the fit. This means fitting the linear part of the astrometric 
orbit (or the O — C diagram) is useless for the parameter determination. 

The astrometric measurements were adopted from "The Washington Double Star Cat- 
alog" (hereafter WDS0), which incorporates a huge database of astrometric observations, 
but only a small one about the properties of the individual components in these systems. 
There was a problem with identifying the eclipsing binaries in them. Altogether more 
than 13800 systems in the northern and southern sky have been inspected. This was the 
number of systems in WDS with 10 or more astrometric observations. From this large 
number of stars only 31 were eclipsing binaries (according to Simbad database). And from 
this 31 stars in most cases there were no or only a little data set in times of minima. 
Systems with larger times-of-minima data sets which have been analyzed are presented 
below in this chapter. In the next section is the brief survey of other systems. 

In a few cases the astrometry and the behaviour of times of minimum light were 
studied, but these two approaches were usually analyzed separately. Such systems are 
for example 44 Boo, QZ Car, SZ Cam, or GT Mus (besides the systems mentioned and 
analyzed below). The coverage of the astrometric orbit is very poor for some of them. 
For SZ Cam only a few usable astrometric observations were obtained, but the LITE is 

"'^http : //ad. usno.navy.mil/wds/ 
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well- defined and also the third hght in the hght-curve solution was detected ( jLorenz et al. 



19981 ). QZ Car is a more complicated, probably quintuple system - the bright component of 
the visual binary consists of two eclipsing pairs (P = 20.7 d and 6.0 d). There are also only 
a few usable astrometric measurements. Also GT Mus is a quadruple system, consisting 
of an eclipsing and RS CVn component. In many other cases, only measurements of the 
times of minima are available, without astrometry. For some others, astrometry without 
photometry, is only available. Other systems where a stromet r ic obs ervations were obtained 
and the LITE is observable or expected are listed in iMayer fl2004 ). 

The only paper on combining the two dif f erent approaches (LITE and astrometry) 
into one joint solution is that by iRibas et al.l (120021 ) . where a similar method (but not 
the same) as described in this thesis was applied to the system R CMa, but where only 
a small arc of the astrometric orbit was available. Besides the astrometry and LITE also 
the pro per motion o r i the long orbit was analyzed. On the other hand one has to note, 
that in iRibas et al.l (|2002| ) the complete astrometric parameters (with proper motions, 
parallax, etc.) were used, while in this thesis only the relative astrometry of the distant 
body relative to the eclipsing pair was analyzed. From this point of view such an approach 
to the combination of LITE and astrometry is unique and has never been published before. 
Generally, such a combined approach is potentially very powerful, especially in upcoming 
astrometric and photometric space missions. 

Three of the systems presented below in this chapter (nara ely VW Cep, C Phe and 
HT Vir) were selected for publication, see IZasche fc Wdl (|2007h (in print) 



4.1 QS Aql 



The first investigated system is QS Aql (KUI 93, HD 185936, HR 7486, HIP 96840). It 
is an Algol-type eclipsing, and also spectroscopic, binary with a period of about 2.5 days. 
Its apparent brightness is ab out 6.0 mag and the spectral type was classified as B5V 
(according to iHolmgrenI 119871 ). 

The star was recognized to be a variable by Dr.Plaskett from The Dominion Astro - 
physical Observatory from photographic plates take n in 1924 a nd 19 25 (jMillmarull928l ). 
The first photometric observations were obtained by iGuthnicki (jl93ll ). Surprisingly, only 
17 times of minima were recorded since then. This is probably due to the relatively high 
brightness of the object, which would saturate most telescopes with CCD detectors (the 
las t one is from Hip parcos) . 



, one IS irom mp parcos ) . 

Guthnick f 193ll ) recognized the eclipsing nature of the star. Some 40 years later, Knipel 



(1l97ll ) discovered a rapid period change, which occurred at about 1964 (his suggestion) 
and was caused by the periastron passage in the wide orbit around the barycenter. The 
period change was so rapid that the eccentricity of the wide orbit must be very high. 
Unfortunately, during the last decade no minimum time was obtained, the last one is 
more than 15 years old. 

The first astrometric observations were secured more than 50 years ago, but their accu- 
racy is questionable. Because both visual components are similarly bright, there could be 
confusion in the identification of the primary and the secondary, and some measurements 
may be shifted for about 180° in the position angle. Especially due to this reason, we have 



Table 4.1: The parameters of QS Aql. 
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neglected all measurements obtained before 1975. Ivlore recent data are more reliable and 
more precise (since 1976 the observations are mostly speckle interferometric). 

Th e only paper considering the astrometry together with the LITE was published by 
IVIayer fl2004 ). The system QS Aql is presented there as one of the systems where LITE 
could be observed together with the astro metric orbit. Also the warning regarding the 
quality of the old data and their 180° ambiguity is given there. 

In Table 14.21 the observed times of minima of QS Aql and the corresponding epochs 
relative to the ephemeris given in Table 14.11 can be found. The algorithm presented in the 
introduction was used to analyze this system combining the astrometry and the times- 
of-minima analysis. The resultant parameters of the distant-body orbit are presented in 
Table 14.11 In Figs. 14.11 and 14.21 the O — C diagram of the times of minima, and the 
astrometric orbit are shown, respectively. The curves in both figures show the model fit 
corresponding to the resultant parameters given in Tabl e I^TTl 

The spectrum of QS Aql was first c lassified as B3 (jivlillmani Il928l ). but already the 
next spectroscopic analysis by (lHilllll93ll ) indicates a later spectral type B5. In this latter 
paper also the first spectroscopic orbit was calculated. The analysis of this SBl-type binary 
results in e = 0.056 ± 0.027, K = (47.31 ± 1.31) km ■ s-\ = (-14.21 ± 0.98) km ■ s-\ 
After then a few spectroscopic investigations of this binary were carried out. The systemic 
velocities were derived: —4.93 km ■ s~^ (lo w confiderice lev el, no spectrosc opic solution, just 
an es timated mean vel ocity of the system, IMillman 1928 ). —13 km • s~^ Lucv fc Sw eenevI 
1971 ). —15.9 km • s~^ ( Batten et al. 19781 ). and the most recent one by Holmgreil (1987.) 




Figure 4.1: An O — C diagram of QS Aql. The individual observations are shown as dots (primary) 
and open circles (secondary) and the curve represents the predicted LITE. All of the measurements are 
photoelectric or CCD ones. 
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Figure 4.2: A relative astrometric orbit of a system QS Aql. Older measurements were neglected. The 
points represent individual observations (black dots), while the solid curve corresponds to the solution 
described in the text and parameters in Table [47T] The straight lines connect individual observations with 
their expected positions on the fitted orbit (red dots). The cross indicates the position of the eclipsing 
binary on the sky, the arrow indicates the direction of the movement on the third-body orbit, and the 
dashed line represents the line of apsides. 



results in (—14.8 ± 0.2) km ■ s^^. As one can see from Fig J4.3l the observed systemic 
velocities are almost constant over the period of the third body, which seems unlikely. But 
it is necessary to take into consideration the error bars (which are not known for some of 
these points) and also the very rapid change in near the periastron and almost constant 
velocity for a decades. 

Using the derived parameters and Eq. I2.15[ also the mass of the third component of 
the system was computed. I f the total mass of the echpsing binary was assumed to be 
Mi2 = 5.9 Mq (according to lHolmgrenlll987l ). one obtains the third mass M3 = 16.5 Mq, 
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Figure 4.3: Systemic velocity variations in QS Aql. The individual points represent derived systemic 
velocities (see details in text), while the solid curve represents the variation on the long orbit described 
by the parameters in Table 14.11 
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large primarily due to the relatively low orbital inclination of the wide orbit. This means 
that the third body is much more massive than the individual components of the eclipsing 
pair, which seems unlikely. One has to take into account the errors of the resultant 
parameters. Due to the relatively large errors of the inclination and mass function, the 
resultant mass could be somewhere between 3.7 and 93.5 M©, which is a very wide range 
of masses. The masses in the lower part of this interval are sufficient to get a reasonable 
lu minosity of the third body. 



Heintze et al.l (119891 ) discussed the spectroscopic observations by iHolmgrenl (119871 ) and 
the light-curve observations and also concluded that the third light is 1.2 times larger than 
the combined light of the eclipsing pair: I3 ^ 1.2 ■ I12, which means that in the bolometric 
magnitude, the third component should be for about 0.2 magnitude brighter than the 
eclipsing binary: Mbois ~ Afboii2 — 0.2 mag. Adopting the spectral types of the primary 
and secondary to be B5V and F3 (i.e. Mboii2 ~ —2.5 mag), the third body should have 
the spectral type B4. If the third star is a main-sequence object, it should have a mass of 
about 5.2 M©. 

This result lies within the range of the masses received from the combined analysis. In 
conclusion, the presented solution is of a low accuracy, mainly due to a very incomplete 
coverage of the astrometric orbit, but leads to an acceptable solution within the limits of 
the errors. 



Table 4.2: The minimum times of QS Aql from photoelectric photometry. 
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4.2 VW Cep 

The eclipsing binary VW Cep (HD 197433, BD +75 752, HIP 101750) was classified as 
W UMa system and in fact it is one of the most often observed and analyzed sys tem. I ts 



magnitude is a bout 7.3 in V filter , but its spectrum is problematic to cl assify. iPoppei 
(11948 ) and also iKaszas et al.l (119981 ) classified the sys tem as K1+G 5 , whil e iPribuUa et al 



(I2OOOI ) proposed the spectral types G5V_+ GSV, but iKaszas et al.l (119981 ) noted that the 



spectral type G5 is inapplicable and iHilll (119891 ) presented the spectral type KOV. Both 
components are chromosphericaly active. VW Cep is rather atypical, because during the 
primary (the deeper one) eclipse the less massive star (the hotter one) is occulted by the 
larger companion (the more massive and the cooler one). 

The first observations of its light variations were done by Schiltl ( 1926 ). Since 1946, a 
large amount of photoelectric observations was obtained. However, the observed minima 
times did not fit the ephemeris due to the LITE and the mass tr ansfer between comp o nents . 
There were many light-time effect studies of this system and iHerczeg &: Schmidtl (1l960l ) 
proposed the presence of a third body with an orbital period of 29 years and an angular 
distance of the third component between 0.5"and 1.2". 

In 1974, the first successful visual observation of the third component was obtained 
and since then, there were 16 observations of it. Regrettably, the observations near the 
periastron passage are missing (the gap in data is from 1991 to 1999). The last two 
measurements {Q = 231.4°, p = 0.702" and 6 = 232.6°, p = 0.695") were obtained and 
kindly sent by Elliot Horch by a speckle camera in April 2007 (priv.comm.). 

The most corn p lete s et of times of minima is in the most recent period study of VW Cep 
by iPribuUa et al.l (120001 ). The first times of minima are from the 1920's and altogether 
1907 minima were collected. From this set of times of minima 313 measurements were 
neglected due to their large scatter (mostly the visual ones). This new minimum-time 




Figure 4.4: An O ~ C diagram of VW Cep using Solution I. The description is the same as in Fig. 14.11 
The bigger symbols are CCD and photoelectric, while the smaller ones are visual. Most of the recent 
visual observations were neglected. The blue dashed line represents the quadratic term and the red solid 
line the quadratic plus the LITE caused by the third body. 




Figure 4.5: An O — C diagram of VW Cep after tlie subtraction of tlic quadratic term, using Solution I. 
The description is the same as in the previous O — C figures, and the red soUd line represents the LITE 
caused by the third component in the system. 



an alysis is based on a larger data set (about 750 times of minima more than were used 
by iPribuUa et al.l ) . see Fig J4.4[ Two new CCD observations of minimum light of VW Cep 
were obtained at Ondfejov observatory. 

The short-term variati o ns wi th the period of about two years 



Kwed (119661) 



(^see e.g 

and iHendry fc Mochnackil (|2000| )) are probably caused by the surface activity cycles on 
the primary component. Due to this activity an unique interpretation of the behaviour 
of period changes is still missing. IPribuUa et al.l proposed a mass transfer (the quadratic 
term) plus the third and the fourth body in the system (two periodic terms). Nevertheless, 
they were not able to explain the O — C diagram in detail. 

Another approach was chosen in this thesis. Especially due to only a few astrometric 
observations (16 measurements from 1974 to 2007) it was decided to explain only the most 
significant effects in the O — C diagram. There were two different approaches used. In 
Solution I. only the third-body orbit besides the mass transfer (the quadratic term) was 
considered, while in Solution II. instead of mass transfer the fourth body on its very long 
orbit and the third body was considered. This approach was chosen especially because 
of the systemic-velocity variations, see below. The astrometric variation with a period of 
about 30 years has been identified with the O — C variation with the same period. 



4.2.1 Solution I. 

The analysis of the times of minima together with the astrometry led to the parameters 
shown in Table 14.31 and the O — C diagram in Fig. 14. 4[ The times of minima, together 
with the curve which represents the LITE and the mass transfer, are shown in this figure. 
After subtraction of the quadratic term, one gets Fig. 14. 5[ where only the LITE caused 
by the third component is displayed. The quadratic term quotient q = —0.756 ■ 10~^° day 
leads to the period change of about 1.98 ■ 10~^day/yr a nd the rate of raass t ransfer from 
the primary component of about 1.30- 10~^MQ/yr (while iPribuUa et al.l (120001 ) derived the 



Table 4.3: The final results: VW Cep, Solution I. and II. The table is divided into three 
parts, in the first one are eleven computed parameters, in the second one the values from 
the literature and in the last one the quantities computed from the previous parts. The 
values of parallax and distance were adopted from the Hipparcos measurements. 



Parameter 


Unit 


VW Cep - Solution I. 


VW Cep - Solution II. 


JDo 


[HJD] 


2437001.5289 ±0.0034 


2437001.4362 ± 0.0025 


P 


[day] 


0.278316241 ± 0.00000011 


0.278315234 ±0.00000012 


Q 


[day] 


(0.756 ±0.032) ■ 10-1° 


- 


P3 


[yr] 


30.04 ±0.47 


29.99 ±0.34 


To 


[HJD] 


2450402 ± 91 


2450366 ± 52 




[deg] 


235.58 ±3.01 


242.53 ±2.89 


e 




0.628 ±0.035 


0.610 ±0.014 


A 


[day] 


0.0117 ±0.0009 


0.0119 ±0.0010 


a 


[mas] 


447.4 ±24.3 


451.4 ±28.1 


i 


[deg] 


30.2 ±4.2 


28.0 ±3.1 


n 


[deg] 


200.0 ±7.1 


192.1 ±5.7 


Ml2 


[Mo] 


1.37 


1.37 


References 




Kaszas et al. (1998) 


Kaszas et al. fl998) 


TT 


[mas] 


36.16 ±0.97 


36.16 ±0.97 


D 


[pc] 


27.7 ±0.7 


27.7 ±0.7 


ai2 


[AU] 


4.30 ±0.41 


4.59 ±0.45 




[Mq] 


0.0112 ±0.0078 


0.0111 ±0.0068 




[Mq] 


0.73 ±0.32 


0.80 ±0.30 



value 1.38 ■ lO-^Mg/yr). 

The fit is not very satisfactory because of the presence of the chromospheric activ- 
ity of the individua l components, or due to the putative fourth component (see e.g. 
PribuUa et alll2000f ). It is evident, that the recent times-of- minima observations devi- 



ates from the predicted fit. This could be caused by a period jump near 1995. In Fig. 14.61 
the astrometric orbit of the binary with the individual measurements and their theoretical 
positions is shown. Regrettably, no observations near the periastron passage are available. 
The curve represents the theoretical orbit according to the parameters given in Table 14.31 
in agreement with the LITE analysis. Also the orbit according to the Solution II. is shown, 
see below. 

The parameters describing the LITE and astrometric variation are in Table 14.31 and 



could be compared to the parameters derived during the previous analysis by lPribuUa et al. 



( |2000| ). Their values for the third-body orbit are: ps = 31.4 yr, e = 0.77, iv = 183°, 
and atotai = 12.53 AU. Our values are in Table 14.31 except for cttotai — 

12.35 AU, and 

as we can see they differ significantly in several parameters. This is due to completely 
different approach describing the O — C variations. Only the period and the amplitude 
of such variation are comparable, but these are th e most importa nt for our combined 
solution. One has also to disagree with the result by iPribuUa et al.l . that the astrometric 
orbit could not be identified with the LITE3 variation from the O — C diagram. As one 
can see, our new results are in agreement with each other without any problems. 

Also the astrometric orbit could be compared with the previously published one. Most 
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Figure 4.6: Relative orbit of VW Cep on a plane of the sky, for a detailed description see Fig. 14.21 Both 
Solution I. and //. are plotted, the black one for Solution I. and the red one for the Solution II.. 



recently iDocobo fc LingI (120051 ) published the following parameters of the astrometric or- 
bit: p3 = 31.0 yr, a = 485 mas, i = 39.3°, and e = 0.68. If one compares these values 
with the new ones (see Table IT3i) . one can see that the differences are slightly beyond the 
limits of errors. 
If the total mass of the eclipsing binary M12 



1.37 Mr;^ was taken fr om (iKaszas et al. 



19981 ) and the parallax tt = 36.16 mas (from iPerryman fc ESAl 119971 ). the distance to 
the system should be only about 27.66 pc, which results in the third-body mass of 
M3 = 0.73 M0. The distant component is about 2.2 magnitudes fainter than the 
VW Cep itself, so its luminosity and also mass should be much smaller than the mass 
of the eclipsing components. Total bolometric magnitude of VW Cep is about 4.7 mag, 
so the magnitude of the third component is about 6.9 mag, which leads to the spectral 
type of about K 3. The typical mass of this spectral type is about 0.75 M0(according to 
Harmaned 119881 ) . which is in an excellent agreement with our result and within its error 
limits. 

Different systemic ve locities v.y wer e found at different epo chs. These values are: v. 



-35.4 ± 1 0) km • s" ^ flPoppeii [l948h . (-f9.8 ±7 ) km • s"^ fiBinnendiikI Il966h. (-8 ± 



1) km-s"^ flHilll 119891 ). and (-16.4 ± 1) km ■ s"^ f lKaszas et al.lll998f ). In the time plot 



[see Fig. 14. 7p one can see the curve which represents the theoretical variation of caused 
by the orbital motion around the common barycentre. Except for the first one data point 
( lPopperlll948l ). the amplitude of the LITE should be much larger (circa 3 times) than was 
computed. Keeping the astrometric amplitude at the value from the fit, this could only 
be achieved by decreasing the inclination to smaller values and by modifying slightly also 
some other parameters. This speculation could only be verified after more accurate and 
larger data set is available. But considering the spectral analysis and efficiency of the RV 
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Figure 4.7: Systemic velocity variations in VW Cep, using Solution I. The individual points represent 
computed systemic velocities (see details in text). 

investigations, this result is not very satisfactory, which is the reason why the different 
approach was also used, see Solution II. below, which describes better the systemic velocity 
variations. 

The star was also measured by the Hipparcos satellite. During its 3-yrs mission there 
were altogether 68 observations obtained, see Fig. 14.81 These are so-called abscissa mea- 
surements and are one- dimensional. It means that only a time of passage through a certain 
main circle was measured, but one cannot derive exactly where on this circle the star really 
was. The position of these circles (measurements) are represented by the small abscissae 
in Fig. 14.81 The observations are connected by a dotted lines with the theoretical positions 
on the sky marked as big points. The theoretical orbit was constructed according to the 
parameters from Table 14.31 As one can see, regrettably only a small part of the orbit was 
measured, so the Hipparcos abscissae measurements are not very useful at all. 

Another task was to derive the parallax of VW Cep using this combined approach. 
Leaving the parallax as another free parameter, one is able to calculate it from the com- 
parison of the angular and absolute semimajor axis (see section [^^ . Using this method, 
most of the relevant parameters remained nearly the same as above, only the inclination 
changed a bit, being about 5° lower. This led to a higher third mass of M3 = 0.99 Mq. 
The main difference was in the parallax, which decreased from 36.16 mas (Hipparcos) to 
33.63 mas. The parallax would shift the distance to 29.74 pc. The new value is only 
about 2 pc higher then the value derived from the Hipparcos measurements. With more 
precise data points and better coverage of the orbit also this result would be better. For 
a comparison with the previously found parallaxes, see the next section. 

4.2.2 Solution II. 

From the analysis of times of minima and using the long period perturbation by another 
distant component instead of the quadratic term, one gets Fig. 14.91 and after the sub- 
traction of the fourth component Fig. 14.101 It means during the computation process 
altogether 14 parameters {A, p^, i, e, u, Q, To, JDq, P, A4, p^, 64, U4, Tq^^) were derived min- 
imizing the xlomb value. 

Applying this approach to the same data one gets about 12 % better result (comparing 
the sum of square residuals). The numerical values for the individual parameters are 
approximately the same (see Table IT3i) . The parameters of the fourth-body orbit are in 
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Figure 4.8: The Hipparcos measurements of VW Cep, see the text for details. 



Table where Mi^rnm denotes for the minimal mass of the fourth body {i^ = 90°). It 
is obvious that the period is about as long as our data set. One can judge that this 
numerical solution is only an edge-on effect, which fits better the most recent data points. 
As one can see, the times of minima since 1995 deviate from the theoretical prediction 
and also an additional period jump should be implemented into the model to describe the 
data points in detail. Nevertheless, this combined approach was chosen because of the 
RV data (see below). In next few years the behaviour of the times of minima will decide 
which solution is the right one. Until that time this approach is just a hypothesis without 
any proof, only for a better description of the systemic-velocity variations. 

The parameters of the third-body orbit according to Solution II. are close to the values 
from Solution I., that the theoretical orbit of the binary on the plane of the sky is almost 
the same, see Fig J4.6[ 

The only effects which are significantly different are the gamma- velocity variations. It 
is shown in Fig. I4.1H where the dashed and the solid line represent the LITE4 and LITE3 
+ LITE4 variations, respectivel y. As one ca n see, this approach gives a much better fit. 
Except for the first data point (|Popperlll948l ) the systemic velocities follow the lo ng-term 
variation and are almost within its errors near the theoretical values. The value by lPopper 



Table 4.4: The parameters of the fourth-body orbit, VW Cep Solution II. 



Parameter 


P4 




Ll!4 


64 




/(M4) 




Unit 


[yr] 


[HJD] 


[deg] 




[day] 


[Mo] 


[Mo] 


Value 


77.32 


2397303 


282.3 


0.561 


0.096 


0.795 


2.64 


Error 


±0.04 


±14 


±2.3 


±0.008 


±0.012 


±0.055 


±0.45 




Figure 4.9: An O — C diagram of VW Cep using Solution II. The description is the same as in Fig. 14.41 
The blue dashed hne represents the LITE caused by the fourth distant body and the red sohd hne the 
final fit LITE3 + LITE4. 



is affected by relatively lar ge error. Th e scatter of the individual RV data points by lPopper 
is larger than those from iBinnendijkl . which could be caused by the combination of two 
different data sets from different instruments and obtained after niore t han 600 orbital 
revolutions (which could shift the ephemeris). iPribuUa fc Rucinskil (120061 ) suggested that 
the scatter of the systemic velocity data points is instrumental, which seems unlikely for 
such a large amplitude. For the final confirmation of w-^ variations a more accurate and 
larger data set is necessary. 

One could also speculate about the possible visual detection of the fourth component. 
This suggested body is bright enough to be visible and its predicted angular separation 
from the system is about 1.1". On the other hand one has to take into consideration that 




Figure 4.10: An O — C diagram of VW Cep using Solution II. after the subtraction of the LITE caused 
by the distant fourth body. The description is the same as in Fig. 14.51 
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Figure 4.11: Systemic velocity variations in VW Cep, using Solution 11. See details in text. 



the period of such a body is not completely covered by data and could be even higher, as 
well as the amplitude could be much higher (and also the angular separation). Nowadays 
there is no potential star for this, only the star BD-l-74 889 shares common proper motion 
and radial velocity, but it is one degree distant. 

Another task was the distance determination. Due to only slight difference between 
the parameters of the third-body orbit from Solution I. and //., also the parallax and 
distance will be approximately the same. The parallax decreased from (36.16 ± 0.97) mas 
{Hipparcos) to (35.85 ± 0.37) mas. This parallax would shift the distance from (27.7±0.7) 
pc (Hipparcos) to (27.90 ± .29) pc. Besides the Hipparcos value, the most precise parallax 
was derived by iHeintzl (119931 ) from trigonometry, resulting in (38.2 ± 1.9) mas. As one can 
see, the values of the parallax determined by Solutions I and II are more precise than any 
of the previously derived par allaxes. For the summary of the previously derived values see 
Hendrv &: Mochnackil fl2000h . 

To conclude, the predicted third body is spectral type K3 with the mass around 
0.73 M0 (according to the Solution I.), or spectral type K2 with the mass around 
0.80 M0 (applying the Solution //.). It is clear, however, that a more complicated model 
will be needed to describe the observed changes completely. Also new times-of- minima 
observations would be helpful to identify the variations in O — C diagram in detail, be- 
cause neither the Solution I., nor the Solution II. are able to describe the behaviour of the 
recent minima observations. This could be described only applying the hypothesis of an 
abrupt period jump. Precise RV investigation (till 2010) would solve the question about 
the nature of the variations in gamma velocity. 



4.3 C Phe 

The system C, Phe is the brightest eclipsing binary with two components of early spectral 
types, exhibiting total and annular eclipses. This is the only eclipsing binary with an 
eccentric orbit included in this study. C Phe (HD 6882, HR 338, HIP 5348) is an Algol- 
type eclipsing binary. Apparent brightness of the system is about 4 .0 mag in V " filte r and 
the spectral types were determined as B6V + B8V (according to Andersenl ( 1983 ). see 



Lingll2004h . It is 



also a comment on the spectral types and Hipparcos measurements in 
a visual triple and SB2 spectroscopic binary. The depth of the primary minimum of the 
eclipsing pair is about 0.5 mag, the period of about 1.7 day. 
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Figure 4.12: Relative orbit of ( Phe on a plane of the sky, for a detailed description see Fig. 14.21 Two 
measurements (the open circles) were neglected. 



It is the visual triple system, while the brightest component is the EB, the most 
distant component is the faintest (some 6"away and with a magnitude of about 8, this 
star is probably not gravitationally bound with the system). The third component is a 
T^'^-magnitude star at a distance of about 600 mas. This is the astrometric component 
and this star is supposed to cause also the LITE variation. 

The first astrometric observation of the third component came from 1930 's and till now 
there were collected 14 observations, but two of them w 'ere neglecte d (see Fig J4.12p . 

Th e unfiltered light curve was observed i n 1950's by Hogg f l95lh . after then by Dachs 



(119711 ) in UBV filters, and the best one by I Clausen et al.l (119761 ) in ubvy filters. In this 
latter paper all relevant parameters of the eclipsing system were derived and also the third 
light was computed. Its value changes from 3%{u) to 8%{y) and the di stant cornponent 
was classified as a spectral type A7 star (the same result was derived by lAndersenI (119831 ) 
on the basis of hi s spect roscopic observations). 



Clausen et al.l (1l976l ) also collected the times of minima obtained before 1975. They 
concluded that n o significant ap s idal m otion is observed. The first apsidal-motion study 



was published by lGimenez et al.l (119861 ). With an updated list of the times of minima one 
is able to conclude that the apsidal motion is definitely presented. It is clearly visible in 
the O — C diagrams shown in Figs. 14.131 and 14.141 Altog e ther 3 6 times of minima use d 
here came f rom t he paper cited above and from Mallania ( 1981 ). Gimenez et al.l ( 19861 ). 
Kviz et al.l (119991 ). The most recent ones are taken from lZasche fc Wolj (l2007h (in print). 
C Phe has one of th e shortest apsidal motions among the eclipsing binaries (see e.g. 



Claret fc Gimeneall993l ) . Due to a low eccentricity, the amplitude of the ef fect is small. For 



an accurate calculation of the apsidal motion rate the method described by lGimenez fc Garcia-Pelayo 



Table 4.5: The parameters of C Phe. 



Parameter 




P 


Vz 


A 


To 


ijj 


e 


I 


{1 






Unit 


[day] 


[day] 


[yr] 


[day] 


[day] 


[dcg] 




[dcg] 


[dcg] 


[Mo] 


[Mo] 


Value 


2441643.7382 


1.6697772 


220.9 


0.0808 


2419900 


97.1 


0.366 


64.4 


33.5 


0.056 


1.73 


Error 


0.0008 


0.0000013 


3.-5 


0.0080 


2500 


2.2 


0.082 


3.0 


4.9 


0.017 


0.26 



( 1l983l ) was routinely used. The eccentricity of the orbit in the echpsing binary is e' = 
0.0107 ± 0.0020, the longitude of periastron ujq = 12.96° ± 5.96°, and the apsidal mo- 
tion rate uj = (0.028 ± 0.001)°/cycle = (6.16 ± 0.20)°/yr, i.e. the apsidal motion period 
U = 58.5 yr. The most recent apsidal-motion analysis is more than 20 years old, made 
by Gimenez et al.l ( 1986 ). but with no LITE and with a smaller set of times of minima. 
The eccentricity by Gimenez was almost the same, but the apsidal motion rate uj was 
0.0373 °/cycle and the angle uq = 13°. 

The approach presented here was a combination of the two different effects. The be- 
haviour in O — C diagram was described as a sum of apsidal motion and LITE contribution 
(O — C) = (O — C)apsid + {O — C)litE) distinguishing the primary and secondary min- 
ima. It means the least-squares algorithm was minimizing the Xcomfe with respect to 12 
parameters in total (v4,p3, i, e, uj, VL, Tq, JDq, P, uj, ujq, e'). 

The astrometric solution based on the combined approach is satisfactory, while the 
older measurements have larger scatter then the recent ones (the old ones are visual and 
the modern speckle-interferometric). Two measurements were neglected, because of their 
large scatter (see Fig l4.12l) . The solution led to the parameters listed in Ta ble 14.51 One 
could compare this new orbit with the previously found one, the most recently iLingI (j2004j ) 
reported the parameters: ps = 210.4 yr, e = 0.348, a = 804 mas, i = 61.9°, 
fl = 33.5°, Go's = 271.7°. It is evident that the new parameters are in very good agreement 
with these ones. The new values imply the mass function of the distant body /(M3) = 
0.056 Mq and with the masses of primar y and secondar y component of the eclipsing 
binary Mi = 3.93 Mq and M2 = 2.55 ]VI0 (lAndersenI Il983l ) . the mass of the astrometric 
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Figure 4.13: The O — C diagram of ( Phe. The apsidal motion curve (the blue one for primary and the 
red one for secondary) is plotted around the (yellow) LITE curve. 
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Figure 4.14: The O — C diagram of ( Phe after subtraction of the LITE. Only apsidal motion curve is 
shown (the blue one for primary and the red one for secondary). 



third body was derived M3 = 1.73 M©. This value corresponds to a spectral type around 
A7, which is in excell ent agreement with the photometric analyses by I Clausen et al.l (119761 ) 
and lAndersenI (119831 ). which result in A7. 



There were also 2 RV investigations by iPopperl (Il970l ) and lAndersenI (119831) . but w ith 
only 2 values of the velocity one cannot do any reliable analysis. In lAndersenI ( 1983 ) is 
also presented that the lines of the third component are also observable in the spectrum 
of ( Phe, but these lines are hardly separable from the binary lines. 

To conclude, ( Phe shows astrometric as well as LITE variations, which are in agree- 
ment with each other. Regrettably, the period of the third-body orbit was not sufficiently 
covered by the data yet, only about 1/4 of the orbit is covered in both methods. Only fur- 
ther precise measurements would prove the third-body hypothesis with higher certainity. 



4.4 V505 Sgr 



Another EB system with apparent changes of the orbital period is V505 Sgr, where the 
third body has been known for more than a decade. It is an Algol-type eclipsing binary 
with a period of about 1.2 days. V505 Sgr (B P 187949, HR 7571, HIP 97849) was classified 
as A2V-I-G5IV spectral types (according to IChambliss et al.lll993l ). with magnitude of 
about 6.5 in V filter. 



The star was discovered to be an eclipsing binary by iHoffmeisterl ( 1l934bl ). Since then a 



lot o f light-curve meas urements and analyses were carried out (for example iLazaro et al. 



20061 ) . The last one by llbanoglu et al.l (I2OOOI ) indicated that the mass ratio of the binary 
is about 0.5 and the contribution to the total light of the binary by the third component 
is about 2.6% in B and 3.6% in V filter. This analysis also yielded the spectral type of 
the distant component to be roughly F6, which is in good agreement with the previous 
analysis by iTomkinI (Il992h. which results in A7. The spectroscopic nature (SB2) was 
discovered by IPoppen (119491 ). 




Figure 4.15: The O -C diagram of V505 Sgr. 



In 1985 an astrometric component was found by iMcAlister et al.l (119871 ) by speckle 



interferometry. The body was 0.3"away from the echpsing pair and after a few years a 
few measurements (16 till now) was obtained. Nowadays it is evident that the distant 
compo nent is moving on its orbit around the EB pair. At the same time also iTomkin 
(119921 ) found the third component lines in the spectrum of V505 Sgr. 

There are several analyses of i ts apparent orbital period ch anges interpreted as the 



LITE due to the third body (e.g. iRovithis-Livaniou et al.lll99ll ). The only paper which 



compares the astrometry and a perio d analysis of O — C deviations from the constant 
orbital period was published by Mayer (fl997i) . Despite existing astrometric measurements, 
there were no attempts to combine these two methods together. The results from different 
approaches were just compared to each other. The main reason why such a combined 
solution is missing, are the differences in parameters, which result from separate solutions 
(see below). 

The O — C diagram is in Fig. 14.151 and the astrometric orbit in Fig. 14.161 The last 
one astrometric measurement was obtained on the 18*^^ October 2005, using 3.6-meter 
CFHT on Hawaii Islands, resulting in p = 0.183(4)", 6 = 218(2)° (kindly sent by Theodor 
Pribulla). As one can see, this point does not fit the theoretical orbit well, but one cannot 
ignore this data point, because it is the only measurement during the last decade and it 
is as precise as the previous ones (the position was obtained after averaging 5 frames). 

As one can see from Fig. 14.151 also the last times of minima in the O — C diagram do 
not follow the theoretical curve, and it is really necessary to observe at least one precise 
minimum of V505 Sgr (the last one is taken from the VSNET database and is not very 



Table 4.6: The parameters from the combined solution of V505 Sgr. 



Parameter 




P 


P-i 


A 


To 




e 


i 


n 


f{M-i) 




Unit 


[day] 


[day] 


[yr] 


[day] 


[day] 


[dog] 




[deg] 


[dog] 


[Mo] 


[Mo] 


Value 


2443750.6866 


1.1828688 


36.86 


0.0085 


2451197 


184.1 


0.802 


195.6 


22.3 


0.0109 


2.76 


Error 


0.0004 


0.0000002 


0.09 


0.0005 


19 


3.2 


0.008 


2.4 


1.9 


0.0024 


0.98 
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Figure 4.16: Relative orbit of V505Sgr on a plane of the sky, for a detailed description see Fig. 



precise, derived only from 9 points). 

The star was also measured by the Hipparcos satellite. Altogether 50 measurements 
were obtained (see Fig. 14.171) . Regrettably, at that time the star was not near its perias- 
tron, so only a small arc of the orbit is covered. It is similar as in the case of VW Cep, 
also the description of the figure is the same. 

The diagrams were plotted according to the parameters from the combined solution 
introduced in Table 14. 6[ It is obvious that the fit to the individual data points is not 
very satisfactory. This is due to inconsistency of the two separate solutions. Only LITE 
solution leads to the 41-years orbit, while the astrometric to 33yr. The angle u differs 
about 60° and the amplitude of the astrometric variation is about 2 times larger than 
one would expect from the LITE analysis. These are the principle reasons why there is a 
doubt of identifying the astrometric and LITE variation to be caused by the same body. 

Because the third body is also visible in the spectra of V505 Sgr, different radial 



velocities of the third component were measured from 1979 to 1989 (see lTomkinlll992l ). 
These measurements together with the predicted variation based on the parameters from 
Table 14.61 are shown in Fig. 14.181 As one can see, there is some systematic increase in 
the radial velocities, but due to only a small part of the period covered, this is not very 
conclusive result (regrettably all the measurements were obtained near apastron). 

From the combined solution, together with the parallax from the Hipparcos, one is 
able to derive the mass of the third bod y. This results in 2.76 M©, or the spectral type of 
about B8 (according to lHarmanedll988l ) . This result is in contradicti on with the pr evious 
spectral analysis, which indicates a spectral type of about F6 (see e.g. lTomkinlll992l ). with 
its typical mass about 1.3 Mq. The only acceptable explanation could be that the third 
component is also a binary. 

In conclusion, the combined analysis of V505 Sgr leads to results which are not in 
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Figure 4.17: The Hipparcos measurements of V505 Sgr, see Fig |4.8l and text for details 



agreement with previous analyses. The spectr al type a r id th e mass of such a body is 
in contradiction with the spectral analysis by iTomkinI ( 1l992l ). This is especially due 
to the inconsistency of the results from the separate LITE and astrometry. Only one 
astrometric measurement was obtained during the last decade, which is not sufficient for 
the precise determination of the orbit. Also in the O — C analysis new precise times of 
minima are needed to prove the LITE variation. The similar situation also apply with the 
radial-velocity measurements. Obtaining the spectra and the radial velocity of the third 
component would be very helpful. 
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Figure 4.18: The radial velocity variations of the third component in V505 Sgr. The solid line represents 
the variation caused by the third body according to the parameters from Table 14.61 
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Figure 4.19: Relative orbit of HT Vir on a plane of the sky, for a detailed description see Fig |4.2l 



4.5 HT Vir 



One member of the visual binary STF 1781 is the echpsing binary system HT Vir (ADS 9019, 
HD 119931, HIP 67186, BD+05 2794). HT Vir is a contact W UMa system, with a period 
of about 0.4 days and the depths of minima of about 0.4 mag. Both visual components 
have almost equal brightness. The third component of the system is brighter than the 
eclipsing binary HT Vir during its eclipses and fainter than it during its maxima. The 
system is apparently about 7.2 mag bright in V filter and the spectral type was classified 
as F8V f accordi ng tolLu et al.ll200lh. 



According to IWalker fc Chamblisd (119851 ) the distant astrometric component was dis- 
covered by Wilhelm Struve in 1830 at a separation of about 1.4" and position angle 240°. 
Since then, numerous astrometric observations were obtained (altogether 277, from which 
275 were used in our analysis) and the orbit is almost completely covered by the observa- 
tio ns (see Fig. 



Baizd ( 119721 ) suggested that the star might be variable. After then, IWalker fc Chambliss 
( 19851 ) obtained a complete light curve of HT Vir and did the first analysis. It indicated 
that both components of the eclipsing pair are almost identical and in contact. The 
temperatures of both components are about 6000 K and the spectral type is estimated 
as F8V. The s ame (combined) spectral type was derived from the spectral analysis by 



Lu et al.l (l200lh . but with a strong con tribution of the third component, 
of the echpsing pai r is M12 = 2.3 My dO'Angelo et aLlliooeh . 



The total mass 



Lu et al.l (I2OOII ) discovered that the distant component is also a binary. They have 



measured the spectra of HT Vir eclipsing pair, and discovered also the lines from the third 
component in the spectra and their RV variations with a period of about 32.45 days. We 
therefore deal with a quadruple system. 
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Figure 4.20: An O — C diagram of HT Vir. The description is the same as in the previous O — C figures, 
all minimum times are the photoelectric or CCD ones. 



Despite the spectral analysis and a large set of astrometric observations, there were 
only a few times of minima published during the last few decades. The main reason is 
the relatively recent discovery of the photometric variability of HT Vir. The first times 
of minima come from 1979. Since then, there were only 31 observations obtained (see 
Fig. I4.20p . Four new observations were obtained, two of them were observed at Ondfejov 
observatory, one by L.Brat and the last one by R.Dfeveny. One unpublishe d observation 
by M.Zejda was also used and four times of minima by M.Zejda published in lZejdal ( 120041 ) 
we re recalculated, becau s e the heliocentric correction was wrongly computed. 



Walker fc Chamblisd (119851 ) published the first rough estimation of the proposed am- 



plitude of LITE from the parameters of the astrometric orbit. Their value (0.18 day) is 
not too far from the present one (0.13 day). 

The final plot of the relative astrometric orbit of HT Vir is in Fig. 14.191 The results, the 
parameters of the orbit around the common barycenter of the system, are given in Table 
14.71 The values of these parameters (A, pa, z, e, a;, fi, Tq, J-Dq, P) were obtained minimizing 

the Xcomb- I 1 

Th e new elements for the astrometric orbit could be compared to these by iHeintz 



(119861 ) ■ which are the following: = 274.0 yr, e = 0.638, a = 1010 mas, i = 42.7°, 
Q = 176.4°, uj^ = 250.0°. As one can see, the period of the new orbit is a bit shorter, 
but the main differences in these values are the angles u and Q. The same fit to the 
astrometric data could be reached with simultaneously transformed values tus — > cus + 180° 
and f2 — i> f2 + 180°. This only means the interchange of the role of the two components. 
This result therefore indicates the incorrect identification of the variable HT Vir in the 
system in our analysis (the variable was s upposed to be the componen t A) and also in the 
WDS catalogue, see WDS notes |. While iPribuUa Sz Rucinskil (l2006h correctly identified 
the variable HT Vir as a B component and A as a single-lined binary. 

If one adopts these parameters to estimate the mass function of the distant pair (mass 
function of the whole pair, not the individual components), one obtains /(M3) = 0.17 Mq. 



^http : //ad. usno.navy.mil/wds/wdsnewnotes_main.txt 



Table 4.7: The parameters from the combined solution of HT Vir. 
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[day] 
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[dcg] 
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Value 


2452722.5040 


0.4076696 


260.7 


0.1274 


2442832 


250.9 


0.640 


45.4 


180.8 


0.169 


2.10 


Error 


0.0050 


0.0000025 


0.5 


0.0026 


61 


0.7 


0.005 


3.7 


2.6 


0.009 


0.11 



This is quite a high value, dictated by the large amplitude of the LITE. With the total mass 
of the primary and secondary M\2 = 2.3 M r:) one gets the t h ird rn ass of M3 = 2.14 Mq. 



The mass of the distant pair is quite high (ID'Angelo et al.l (120061 ) derived the mass for 



some 50% lower, M3 = 1.15 Mq), but note that also this object is a binary and we do 
not know the individual masses. From the spectroscopic observations (to remind, it is a 
SBl-type binary), one is only able to estimate the mass function of the components, or 
some upper limit for one of them (we do not know the inclination). The present result M3 
is the total mass of the SBl pair M3 1 + M3 2; the limit for the invisible-component mass 
M3 2 ■ sin(z') = 0.075 Mq. If the the coplanar orbit is assumed, high difference in masses 
would arise, one component should be much more luminous and also more luminous than 
the eclipsing pair itself, which is not the case. In fact the whole system is not coplanar 
(see e.g. i = 315.5° and the inclination of the EB close to 90°). If one assumes two 
approximately equal masses, there is a problem with the luminosity, because the distant 
pair has to be roughly as luminous as the eclipsing pair. This could only be satisfied if 
one component is underluminous or degenerate. 

One has to take i nto co nsideration also the comment on the light-curve solution by 



Walker fc ChamblissI (119851 ). Using the Wood's model, they discovered that if the third 
light L3 is fixed to be the equal to the light from the distant visual component (it means 
L3 = 0.5), the solution of the light curve is unrealistic. To conclude, the system could 
be much more complicated than the approach that was used here. There may be some 
additional component (s) or the distant pair is composed from evolved stars, away from 
the main sequence. Especially because of the resultant mass and luminosity of the distant 
pair, the body causing the astrometric variation is probably different from the one causing 
LITE, but this conclusion will be proven only if also the nonlinear part of the O — C 
diagram is covered. 



4.6 The problematic case: V2388 Oph 



Another system where the astrometric orbit is known and also the set of times of minima 
is available is V2388 Oph. The contact eclipsing binary system V2388 Oph (FIN 381, 
HD 163151, HR 6676) is P Lyra e type (according to the Simbad database), or more 
likely W UlVIa type (according to iRodriguez et al.lll998l ). Its orbital period is about 0.8 
days, apparent brightness about 6.3 mag in F filter and the spe ctral type was classified 



as F5Vn (acco r ding to Hipparcos Catalogue, iTuron et al.l Il993l ). or F3V according to 



Rucinski et al.l (120021 ). The depths of its minima are 0.3 and 0.25 mag for primary and 



secondary, respectively. 



The first astrometric observation of the third component was obtained in 1959 (iFinsen 



1963! ) in a distance of about only 100 mas. During the next decades there was recognized 
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Figure 4.21: Relative orbit of V2388 Opli on a plane of the sky, Solution I. Three measurements were 
neglected due to their large scatter. 



ra pid mo v emen t of this component around the primary. A prehminary orbit was calculated 
by iBaisd ( 119881 ) with a period of about 8.3 yr, semimajor axis 0.09"and eccentricity 0.29. 
Nowadays orbit is a little bit different (pa = 8.9 yr, a = 0.09", e = 0.33). It is evident, 
that the movement is very rapid and since its discovery the body revolved 5 times around 
the eclipsing pair. 

The photometric variability was discovered by iRodriguez et al.l (119981 ). but the vari- 
ability is also evident from the Hipparcos observations. The ubvy ligh t curv es were 
obtained and analyzed. With the RV analysis made by iRucinski et al.l (120021 ) one is 
able to get the complete picture of the system. The mass ratio from the photometry 
is 0.27, but more precisely from spectroscopic analysis q = 0.186, the minimum mass 
(Ml + M^ jsin^i = 1.93 Mf7>. According to ubvy analysis the system was classified as 
F5Vn, but lRucinski et ahl on the basis of their spectroscopic observations suggested slightly 
earlier spectral type F3V. The magnitude difference between astrometric components is 
Am = 1.80 mag and the contribution of the third component to total luminosity i s circ a 
20%. Mass of the third component is M3 = 1.36 Mfr^ accordin g to lD'Angelo et al.l (120061 ). 
The EB system is SB2-type (according to lRucinski et al.ll2002l ) also with the third compo- 
nent visible in the spectra. The mean radial velocity of the third component VC^, 3 = 30.64 
km-s~^ significantly differs from the center-of-mass velocity of the binary, Vq = 25.88 
km-s~^, which could be caused by the motion on the 9-yrs orbit. In the same paper was 
mentioned an unexplained cross-talk, which means variation in radial velocities of the dis- 
tant companion in phase with the period of the eclipsing binary. Having only this one RV 
data point one could not determine any of the relevant parameters of the long orbit. 

Rough estimation of the proposed O — C variation and its magnitude was done by 
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Figure 4.22: An O — C diagram of V2388 Oph, Solution I. All minima times are the photoelectric or 
CCD ones (the first 5 of them are from the Hipparcos mission). 



Rucinski et al.l (120021 ) from the astrometric orbit. Their analysis results in A = 0.011 day. 
Since the discovery of the photometric variability of V2388 Oph the minimum light of the 
star was observed only 18 times. But the individual data points show different variation, 
with quite different amplitude, but mainly with very different period. 

The astrometric orbit is plotted in Fig J4.21l and the final fit is satisfactory, with the 
resultant parameters given in Table 14.81 - Solution I. The O — C diagram is in Fig j4.22[ 
As one can see, the fit is unacceptable, but the method applied here was the same as in 
the previous cases. 

Where the problem could be? The crucial point at the first time is to compare the 
results from these different approaches and take it into the consideration. If these two 
results are incompatible (as in this case) this combined approach is unusable. And this is 
the case of V2388 Oph. 
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Figure 4.23: An O — C diagram of V2388 Oph, Solution I. - only the quadratic term was used. 



Table 4.8: The final results: the case of V2388 Oph, Solution I. and 11. The description is 
the same as in Table 14.31 
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Using only — C analysis to the set of times of minima, one gets the period = 5.1 yr, 
e = 0.56 and lo = 102.4°. It is evident that the astrometry leads to the different set of 
parameters and the combined solution should be nonsense. Another solution could be the 
very short one, with the period ps = 2.99 yr, but this is only the hypothesis, because the 
variation in O — C diagram is not covered very well and this is only a sampling frequency 
of the individual data points. 

The astrometric orbit is well defined, but there is a question about the accuracy of the 
individual times-of- minima data points. Without the input data (the rough photometry), 
one can doubt, if all the measurements are accurate enough or some of them could be 
neglected. Interesting is the sequence of 4 times of minima, which are rising up near the 
epoch -1000 (3 primary and 1 secondary). Is this the real effect in O — C diagram, or is 
it just the real scatter of the measurements? The difference is about 0.01 day, or circa 14 
minutes. This is quite large to be only a scatter, but one does not know the conditions 
during the observation, etc. 

If one decide to neglect one data point - the secondary minimum time near the epoch 
0, and include also the quadratic term in the ephemeris, one will get Solution II. - see Fig. 
I4.23[ Using the combined approach also the astrometric orbit could be plotted (see Fig. 
14.241) and as one can see, the difference between Figs J4.2l1 and W?IM is not so significant. 
The main difference is in the — C diagram, where only the quadratic term arises and 
no LITE is presented. This means that the orbit is just face-on, i.e. the inclination is 
very close to 180°. See the resultant parameters of such fit in Table The parameters 
of LITE are not very convincing, because of inclination is almost 180°, but the LITE is 
necessary to compute, because the amplitude of astrometric variations is computed from 




Figure 4.24: Relative orbit of V2388 Oph on a plane of the sky, Solution II. 



the amplitude of LITE. This means the mass of the third component was derived precisely, 
while the mass function of such a body is very inaccurate. 

To conclude, it is difficult to decide which solution is the right one. Only further 
data points, especially times of minima, would confirm the 9-yrs variation in the O — C 
diagram. Due to very short period of the third body the shift in the times of minima 
should be evident after a few months of observations. Also measuring the RV curve and 
the analysis of the systemic velocity could be very helpful, because the period is short and 
the last one was carried out more than 6 years ago. 



4.7 Other systems 



In this section are presented the systems which were found to be EBs as well as members 
of the visual binaries. The limitation about the number of times of minima, which was 
presented in Introduction to Chapter HI does not play a role. Only the systems analyzed in 
detail above were omitted. The systems presented here were found by scanning the objects 
in the WDS catalogue and trying to identify the EBs in this sample of stars. This survey 



i s slig htly following the paper on "Eclipsing binaries in multiple-star systems", IChambliss 
(I1992I ). The number of such s ystems has grown rapidly since then, but the main difference 
is the selection criterion. In IChamblisd (Il992l ) are presented all of the m ultiple systems 



with eclipsing binaries, which were known for the author. IChamblisd mentioned that 
80 EBs are known to be components of the multiple-star systems and 37 of them were 
presented in more detail. For this thesis there were sele cted only these sy stems, which 
were discovered to be EBs and astrometric variables after IChamblisd (jl992l ). the systems 
which have sufficiently large data set in both methods to do the simultaneous analysis, or 



the systems for which the new astrometric orbit was calculated for the first time. Some of 
the presented systems were also included due to their misidentification as EBs. 

4.7.1 HD 123 

HD 123 (V640 Cas, HR 5, STF 3062AB) is an eclipsing binary which spectral type was 
classified as G5V. Its V magnitude is of about 5.93, but there were only a few times of 
minima found in the published literature, no photometric analysis was found. The eclipse 
observations are questionable and recent measurements indicate possible misidentifications 
of the star as eclipsing binary. On the other hand the as trometry coyers w hole orbit. 



Altogether 572 data points were obtained during 170 years. ISoderhjelml (1l999l ) computed 



the orbital parameters, the period about 107 yr and angular semimajor axis about 1.4' 
4.7.2 HD 1082 

HD 1082 (V348 And, A 1256AB, HIP 1233) is an Algol-type EB, which spectral type was 
classified as B9V. Its apparent magnitude is 6.76 in V filter. The same situation as in the 
previous case also apply here, there were neither no times of minima nor the photometric 
analysis found in literature. The astrometric orbit is covered by 61 data points obtained 
during 93 years and co vering the ran ge from 4 to 223 degrees in Q. From these data the 



orbit was calculated by lOlevid (l2002l ). resulting in = 138 yr and a = 150 mas. 



4.7.3 HD 4134 

HD 4134 (V355 And, STF 52AB, HIP 3454) is also an Algol-type EB with spectral type 
classified as F5 and the magnitude V = 7.69 mag. No times of minima were obtained. 
Astrometry covers only 20° with 51 measurements observed in the last 170 years, the orbit 
was not computed. 

4.7.4 HD 10543 

HD 10543 (V773 Cas, BU 870AB, HR 499) is an Algol-type EB with spectral type A3V 
and the magnitude V = 6.21 mag. The orbital period of the eclipsing pair is about 1.3 
days, but only one time of minimum was observed. The astrometry covers about 80° 
with 79 obse rvations made durin g 120 years. The astrometric orbital parameters were 



calculated by lPopovic fc Pavlovid (119951 ). resulting in period about 304 yr and semimajor 
axis about 1". 



4.7.5 HD 12180 

HD 12180 (AA Get, ADS 1581 A, HIP 9258) is W UMa type EB, sp F2V, V = 7.22 mag, 
and orbital period about 0.54 d. There were more than 200 times of minima obtained 
during the last 40 years, but with no significant LITE variation. Also the astrometric 
observations, which were obtained during more than 200 years, do not show any evident 
variation and any orbital solution could be found from this data set. 



4.7.6 HD 14817 

HD 14817 (V559 Cas, STF 257AB, HIP 11318) is one component of the visual binary 
STF 257AB. It is the eclipsing binary of Algol-type, as well as spectroscopic binary, spec- 
trum classified as B8V, apparent brightness of about 7.02 mag in Y filter and orbital 
period of about 1.58 day. There were 7 times of minima observed since 1971 to 1991. Due 



to its very long orbital period, about 836 yrs (see e.g. iHartkopf et al.ll200ll ). only about 
one third of the orbit is covered by the observations (101 observations and the change in Q 
is about 100°). The astrometric measurements are available since 1830 and the periastron 
passage occurred in 1932, so the part of the orbit near periastron is sufficiently covered. 
Regrettably, in that time the minima times are missing. 



4.7.7 HD 18925 



HD 18925 (7 Per, 23 Per, HJ 2170A, HR 915) is an Algol-type EB with spectral type of 
about G8III and the magnitude V = 2.95 mag. Astrometric observations (alto gether 67) 
were obtained during 65 years and the orbit was calculated. iPourbaix (l2000f ) published 
the parameters, = 14.6 yr and a = 144 mas, and the inclination i = 90.6°, so the orbit 
is just edge-on and the astrometric observations are only "in the line" . The position of the 
orbit indicates that the occultations and ecli pses may happen. T hese were predicted and 



successfully observed in September 1990 (see iGriffin et al.l (1l994l ) for details). The system 



is similar to (3 Aur (see below). Therefore, this object is not suitable for the simultaneous 
analysis. 



4.7.8 HD 19356 

HD 19356 (Algol, p Per, LAB 2Aa, HR 936) is well-known prototype of the Algol-type 
binaries. Its spectral type is B8V and apparent brightness 2.12 mag in V. The time of 
minimum brightness was first measured by Montanari on 8 November 1670 (although 
known from historical times). Nowadays set of times of minima is really large, about 1400 
observations, covers a few centuries, but the detailed description of the O — C diagram 
is still missing. The system is rather complicated, but the distant component with the 
orbital period about 1.8 yr discovered firstly on the basis of the radial velocity variations 
was found in 1973 by a speckle camera and the o rbit of this com ponent is now well 



established (a = 94.6 mas and e = 0.23, according to lPan et al.lll993l ) 



4.7.9 HD 24071 

HD 24071 (DUN 16, HR 1189, HIP 17797) is probably p Lyrae type star, its apparent 
brightness is about 4.2 mag in V filter and spectrum classified as B9V. The star is hardly 
measurable, because there are together 4 stars very close each other (only 5"distant) and 
it is not clear, if all these components belong to the system. The astrometry was obtained 
in 1826 for the first time, there were 80 measurements in total, which reveals the change 
in position angle of about 15°. The orbit computed according to these data is not very 
conclusive (the orbital period more than 5000 yr). 



4.7.10 HD 25833 



HD 25833 (AG Per, STT 71AB, HIP 19201) is an Algol-type EB spectral type B5Vp and 
the relative magnitude V = 6.69 mag. There were 101 times of minima, collected from 
the published literature. These minima were obtained from 1920's till now. AG Per is 
one of the most typical apsida l-motion systems , which has been analyzed for the apsidal 



motion several times (see e.g. IWolf et al.ll2006l ). One could also apply the hypothesis of 



combining the apsidal motion and the LITE into one joint solution ( similar to the C Phe 



case) . Also the precise light curves were measured and analyzed (see IWoodward fc Koch 
19871 ). The main problem arises with the astrometry. Altogether 38 measurements cover 
more than 30 degrees in 6, but there is no evident periodicity and the orbit could not be 
constructed from this data set. 



4.7.11 HD 29911 

HD 29911 (V592 Per, COU 1524, BD+39 1054) is a /3-Lyrae EB with spectral type clas- 
sified as F2 and its apparent magnitude V = 8.37 mag. There was only one time of 
minimum obtained. The astrometry covers only 19° with 18 data points obtained during 
26 years. The plot is in Fig J4.25[ where the theoretical orbit is also shown. Its period is 
about 117 yr and semimajor axis of about 230 mas, it was computed for the first time. 

4.7.12 HD 36486 

HD 36486 {6 Ori A, 34 Ori A, HEI 42Aa, HR 1852) is an eclipsing binary, sp 09.511, 
V = 2.23 mag, and an orbital period 5.7 days. Only 9 times of minima were found in 
literature, but these minima do not show any significant LITE variation (more probably 
apsidal motion). On the other hand, there is significant motion on the plane of the 
sky, while the astrometry was first obtained in 1978 and since then 38 observations were 
obtained (see Fig J4.26p . The orbit is only a preliminary one and has not been published 
till yet. The period of the orbit is about 313 yr and the semimajor axis 280 mas. 
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Figure 4.25: Relative orbit of V592 Per on tlie plane of the sky. 
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Figure 4.26: Relative orbit of 5 Ori on the plane of the sky. 



4.7.13 HD 38735 

HD 38735 (V1031 Ori, MCA 22, HR 2001) is an Algol-type detached system, V = 
6.06 mag, sp A4V, period about 3.41 d. There were 9 times of minima found in liter- 
ature. The orbit of the binary is shown in Fig J4.27[ It consists of only 20 observations 
obtained from 1980 to 1997. This orbit was not published yet and is only a preliminary 
one. Its orbital period is about 92 yr and the semima jor axis about 0.18 ". But accord- 



ing to the RV measurements by I Andersen et al.l (Il990l ). the orbit should be much larger. 



and the period about 3700 yr. The third-component lines were observed in the spectra 
of V1031 Ori and radial velocities on the 92 yr orbit would be much larger than mea- 
sured. Because the orbit is covered by data only very poorly, only further astrometric 
observations, as well as precise radial velocity investigation will reveal the nature of the 
system. 



4.7.14 HD 40183 

HD 40183 (/3 Aur, 34 Aur, HR 2088) is an Algol-type EB, apparently bright about 1.9 mag 
in V filter and its spectrum was classified as A2IV. It is one of the brightest and nearest 
spectroscopic as well as eclipsing binaries, but due to its high brightness only a few obser- 
vations were done. Altogether 23 times of minima were measured over the whole century. 
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Figure 4.27: Relative orbit of V1031 Ori on the plane of the sky. 



Photometric (see lJohanserull97ll ) and also spectroscopic (see lNordstrom fc Johansenlll994j ) 



analyses were published. Similarly to the previous case 7 Per, the astrometric orbit could 
be identified with the eclipsing binary orbit. Therefore, this system does not belong to 
this survey, but it is of big importance for the present EB knowledge. Altogether 28 data 
points sufficiently cover the whole 4-day astrometric orbit. The EB components could be 
resolved, because the system is relatively close (about 24 pc). Detailed description of the 
technique used (interferometry with The MARK III long-basel ine optical interferometer 
on Mount Wilson) and the analysis is in iHummel et al.l (119951). The para meters of the 
orbit from interferometric measurements were compared by iHummel et al.l with the pre- 
viously found values from photometry and spectroscopy. The different approaches lead to 
the same results within their respective errors. Also the determination of the distance to 
this unique binary from four independent methods gave the comparable results. 



4.7.15 HD 57061 

HD 57061 (r CMa, 30 CMa, FIN 313Aa, HR 2782) is the brightest star in the open 
cluster NGC 2362. It is a /? Lyrae-type EB, period about 1.28 d. r CMa is also a 
spectroscopic binary with an orbital period of about 154.9 day and the EB is probably the 
main component of the SB. This interesting system therefore contains both the longest 
period spectroscopic binary and the shortest perio d eclipsing binary known among the O - 



type stars. The system was precisely analyzed by Ivan Leeuwen fc van Gendered (119971 ). 
This triple system is one member of the visual binary FIN 313Aa, which has been measured 
32 times since 1951. The change in position angle is only about 15°, so any orbital solution 
is acceptable. 



4.7.16 HD 66094 

HD 66094 (V635 Mon, A 1580AB, BD-08 2186) is an Algol-type EB with primary star 
classified as a spectral type F5 and apparent brightness of about 7.31 mag in V filter. A 
lot of times of minima were collected (altogether 113), but these data points follow the 
linear ephemeris without any indication of the proposed LITE. On the other hand the 
astrometric orbit is defined very precisely. 23 data points measured over the century cover 
about a half of the orbit, and the analysis results in 160- yrs orbit, with the semimajor axis 
of about 280 mas. 



4.7.17 HD 71581 

HD 71581 (VV Pyx, B 2179AB, HR 3335) is an Algol-type EB, spectrum AlV, brightness 
6.58 mag in Y and orbital period about 4.6 days. There were 7 times of minima found in 
literature (1976 - 1983), but these data show very long apsidal motion (in order of decades 
or centuries). The astrometric orbit is also covered only very poorly (11 observations 
obtained during 38 years show the change in position angle of about 13°). 



4.7.18 HD 74956 



HD 74956 (5 Vel, HR 3485, HIP 42913) is an Algol-type eclipsing binary classified as AlV 
spectral type, with V = 1.95 mag. The star was discovered to be a photometrically 



variable in 1997 (see lOtero et al.l (120001 ) for details), the period of such variation is about 
45 days. Altogether 8 times of minima were collected, but these data indicates very 
long apsidal motion (on the timescale of centuries). Astrometric orbit consists of 37 
measurements, which define th e orbit with the period of about 142 years and the semimajor 
axis of about 2" (according to lAlzner fc Argyld 120001 ) . The whole system is in fact more 
complicated, consists of two proper motion pairs (2"and 6") separated by 69". Also the 
primary component was resolved as a double star interferometrically. We therefore deal 
with a quintuple system (at least 5 components). 



4.7.19 AC UMa 

AC UMa (ARC 21B, BD+65 671B, AG+65 453) is an Algol-type EE, spectrum classified 
as A2, brightness 10.3 in V filter. The orbital period is about 6.85 days. There were 69 
times of minima found in literature and there could be some variation in order of decades, 
but this is only hypothesis, larger data set is needed. The astrometry is shown in Fig. 
14.281 only 10 observations during 106 years were obtained. As one can see, only a linear 
part of the orbit is covered by data, so one cannot derive the parameters of the orbit 
precisely. This leads to extremely long period about 1200 yr, which could be even higher. 

4.7.20 HD 82780 

HD 82780 (DI Lyn, A Hya, STF 1369AB, HR 3811) is an Algol-type EB, classified as 
F2V. Its magnitude is = 6.76 mag. There is only one time of minimum measured and 
the small arc of the astrometric orbit was observed during 22 years, covering about 30 
degrees. No acceptable solution could be found. 
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Figure 4.28: Relative orbit of AC UMa on the plane of th.e sky. 



4.7.21 HD 91636 

HD 91636 (TX Leo, 49 Leo, STF 1450AB, HR 4148) is an Algol-type EB. Its spectrum 
was classified as A2V and its apparent brightness is about V = 5.67 mag. There were 6 
times of minima observed since 1930. The astrometric data set is much larger, about 132 
measurements secured during the last 180 years, but the motion is undetectable. 



4.7.22 HD 101205 

HD 101205 (V871 Cen, I 422AB, HIP 56769) is a /? Lyrae type EB with its spectrum 
classified as 08V and the brightness of about V = 6.49 mag. There is a brief paper on 
the photometric o bservations of V871 Cen, together with a minimum time derived (see 



Mayer et al.l 119921 ). The astrometry secured during the last 90 years reveals the change 



in 9 about 20°, but no acceptable solution could be found (these data lead to an orbit of 
period about 4500 yr). 



4.7.23 HD 101379j 

HD 101379j (GT Mus, 12 Mus, B 1705AB) is an eclipsing binary, with spectrum classified 
as G2III and the brightness V = 5.17 mag. The astrometric data were obtained during 
60 years an d cover about 130° of the orbit. This SBl-type spectroscopic binary was 
analyzed by iParsond (120041 ). The orbital period of GT Mus is about 56 days, which is 
different from the period of the SBl binary. Therefore, the component A is the EB, while 
B is the SB. N o minima were derived. The astrometric orbit has a period circa 91 years 
( lParsonsll2004l ). Some observations indicates that one of the components is RS CVn-type 
star. 



4.7.24 HD 103483 

HD 103483 (DN UMa, 65 UMa A, HR 4560) is an Algol-type EB, sp A3Vn, V = 6.54 mag, 
orbital period 1.73 days. There were only twelve times of minima found in published 
literature (the first ones from 1979). The astrometric orbit is covered sufficiently, the 
first astrometric observation came from 1 908, and the parame ters of the orbit are known 
(ps = 136.5 yr, a = 230 mas, according to lAristidi et al.lll999l ). 



4.7.25 HD 110317j 

HD 110317j (VV Crv, STF 1669AB, HIP 61910) is an eclipsing binary with the spectrum 
classified as F5IV and the brightness of about V = 5.27 mag. There were no times of 
minima found in the literature. The astrometric data set consists of 156 measurements 
secured during 180 years, which yielded a change in 6 of about only 14°. 



4.7.26 HD 114529 

HD 114529 (V831 Cen, SEE 170AB, HR 4975) is the p Lyrae system, sp B8V, V = 
4.58 mag, orbital period of about 0.64 d. No published minima were found. The astro- 



metric orbit was derived according to 40 observations secured dur ing the last 100 years, 
resulting in p3 = 27 yr and a = 185 mas (according to lFinsenlll964l ). 



4.7.27 SAO 45318 

SAO 45318 (ET Boo, COU 1760, HIP 73346) is a /5 Lyrae eclipsing binary, spectral type 
F8. Its apparent brightness is V = 9.09 mag. There were found a few times of minima, 
covering the last 5 years. There is possibly some variation in the O — C diagram, but its 
amplitude is only about 0.001 days and the period about 1.25 years. On the other hand the 
astrometric measurements were obtained since 1978 till 19 99, altogether 20 observations 
show the change in 6 about 40°. The orbit was derived by ISeymourl (l200ll ). resulting in 
period about 113 yr and angular semimajor axis 261 mas. 



4.7.28 HD 133640 

HD 133640 (i Boo, 44 Boo, STF 1909AB, HR 5618) is a well-known EB of W UMa type, 
spectral type GOVnv and brightness of about V = 4.76 mag. It is quite a complicated 
system, consists of more than three components. Many times of minima were observed 
during the last 90 years, but the detailed description of the behaiour of these minima is 
still missing (mass transfer + LITE ?). It was found to exhibit flares as well as to be 
an X-ray binary and also many analyses in this part of spectra were obtained. The large 
astrometric data set consists of 753 observations secured during the last 223 years, and 
covers the range of position angle from 2 40 down to 57 deg rees. The orbit has period of 



about 206 yr and semimajor axis 3.8^Ysee ISoderhjelmlll999l ). 



4.7.29 HD 148121 

HD 148121 (V1055 Sco, B 872AB, HIP 80603) is p Lyrae EB with spectral type classifled 
as G3V and brightness V = 8.64 mag. There were no times of minima found in the 
literature. Astrometric measurements were obtained 12 times during the last 70 years 
covering about 15° in position angle. The orbit was not derived. 



4.7.30 HD 157482 



HD 157482 (V819 Her, MCA 47, HR 6469) is an Algol-type EB, spectrum analyzed as 
F9Vn and its magnitude is about 5.57 in V fllter. The EB pair is orbiting around the 
common center of mass with the third component on the 5.5 years orbit with eccentricity 
0.67, LITE is evident. This is the only system where the LI TE was analyzed together with 



the other methods, namely the interferometry and RV (see iMuterspaugh et al.l 120061 ). 



4.7.31 HD 163708 

HD 163708 (V1647 Sgr, HIP 88069) is an Algol-type EB with the spectrum classifled 
as A3III and the relative brightness V = 6.8 mag. A few dozens of times of minima 
are available, showing very slow apsidal motion (in order of centuries). The astrometric 



measurements were obtained 15 times during 170 years and covering about 14 degrees in 
position angle, see Fig |4.29[ The orbit was computed first time and was not published 
yet, but the result is not very convincing due to poor coverage of the orbit by data points. 
The period is about 1200 yr and semiamplitude about 7.7". 



4.7.32 HD 174932 

HD 174932 (COU 510, BD+24 3555, SAO 86519) is a member of visual binary COU 510. 
The system also does not belon g to th is list, because the star was incorrectly classified as 
an eclipsing binary by ICouteaul (1l972l ) and designated as JZ Her. HD 174932 itself is not 
a variable star it was mixed up with the close eccentric eclipsing variable HS Her. The 
wrong comment in the WDS notes will be soon corrected. 



4.7.33 HD 178125 

HD 178125 (18 Aql, Y Aql, HEI 568AB, HR 7248) is a spectroscopic variable, spectral 
type B8III, brightness about 5.07 mag in Y filter. The astrometric orbit of Y Aql leads 
to the period of about 58 yr. In fact, the star does not belong to this survey, because it 
is probably not eclipsing, but rather an ellipsoidal variable (recent observations and also 
data from Hipparcos indicate this possibility). 



4.7.34 HD 184242 

HD 184242 (V2083 Cyg, A 713AB, HIP 96011) is an Algol-type EB, spectral type A3, 
apparent brightness V = 6.88 mag and orbital period about 1.9 days. There were no 
times of minima found in published literature. The astromet r y cov ers about 70° during 
the last century. The orbit was computed by ISeymour et al.l (120021 ) . resulting in period 
about 372 yr and angular semimajor axis about 498 mas. 
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Figure 4.29: Relative orbit of V1647 Sgr on the plane of the sky. 
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4.7.35 HD 195434 

HD 195434 (MR Del, AG 257AB, HIP 101236) is an Algol-type EB, spectral type classified 
as KO and V = 11.01 mag. There were found a few times of minima, covering 2451700- 
2452100 HJD and astrometry from 1902 to 2001, with the change in position angle of 
about 15°. This leads to an orbit with period of about 6500 yr. 



4.7.36 HD 201427 



HD 201427 (BR Ind, HU 1626AB, HIP 104604) is an Algol-type EB, spectral type F8V 
and with its apparent brightness of about 7.1 mag in V filter. Astrometry was obtained 
since 1914 till 2001, when the position angle has changed from 208 down to 124 degrees. 
These data lead to the orbital parameters ps = 167 yr and a = 894 mas (according to 
Seymour et al.ll2002l ). No times of minima were found in literature. 



4.7.37 HD 217675 



HD 217675 (o And, 1 And, BLA 12Aa+WRH 37AB, HR 8762) is a pulsating Be star, as 
well as a shell star. Its apparent brightness is = 3.63 mag and spectrum was classified 
as B6IIIpe. The syste m is more complicated, consisting of at least 4 components (see e.g. 
Pavlovski et al.l 119971 ). while the visual triple con sists of Aa-B co r apone nts. Both orbits 
(the longer one with period 68.6 vr . according to iHartkopf et al.l (1l996l ) and the shorter 
one with period 8.9 yr, according to 



Olevic &: Jovanovic 



1999f ) were observed and derived. 



Some authors (see e.g. ISchmidtlll959l ) published the light curves with the possible eclipsing 
behavior of the star, but nowadays it is rather improbable for the star to be an eclipsing 
binary. The photometric variability is probably due to the variability of the shell around 
the star. 

The survey of 37 systems with eclipsing components in visual binaries is not a complete 
one. On the other hand it could be taken as a representative sample of the most interesting 
ones, because these systems have the largest data sets in astrometry and some of them 
have also the times-of-minima observations. 

Scanning the WDS catalogue and trying to find the eclipsing binaries in this sample 
there were found a lot of variable stars (according to Simbad catalogue). From this sample 
of variable stars there could be a significant number of eclipsing binaries, but only future 
photometric observations would reveal the nature of this variability. Also a few systems 
mentioned above are classified in Simbad as variable stars or ellipsoidal variables, although 
they are EBs. Some of the binaries were also find to be wrongly identified. 



Chapter 5 

Discussion and conclusions 



The method of period analysis of echpsing binaries and its modifications were presented. 
The method of O — C diagram analysis is not new, but new aspects were also included 
into the code. The possibility that the third body resulting from the LITE analysis is also 
detectable via astrometry was discussed. With this assumption the modified algorithm of 
simultaneous solution of LITE and astrometry was presented. 

The theoretical explanation of the effect and the method used is presented in chapter 
El It deals with the relative astrometry only and also the parallax of the system was 
assumed as a priori known (mainly from Hipparcos satellite). The distance is needed for the 
transformation between the angular and absolute semimajor axis in both methods. Also 
the principal limitations for both methods are presented in chapter O These limitations 
have to be considered, especially when one has only very poor data in one of the methods. 

LITE systems: In chapter [3. II there was presented the application of the LITE analysis 
on eleven particular systems. These systems have never been studied for the presence of a 
third component and LITE hypothesis is able to describe their long-term minimum times 
behavior. On the other hand there were neither detailed spectroscopic, nor photometric 
analyses of these systems and the third body hypothesis presented here cannot be proven. 

Although the number of systems, where the astrometric orbit together with LITE is 
known, is growing steadily, in most cases only very limited coverage of the orbit, both in 
astrometry and times of minima is available. Especially due to this reason the combined 
analysis of these systems is still difficult. There were found only a few appropriate can- 
didates for such an analysis. These cases were studied and discussed in detail and the 
principal limitations of the method were pointed out. On the other hand the method itself 
is very powerful and efficient. It could be even modified for the estimation of the distance 
to the suitable kind of binaries. 

During the last decade a few papers combining the approach of simultaneous solution 
of radial velocities, spectral analysis, astrometry, Hipparcos measurements or LITE were 
published. Besides the systems mentioned in the introduction (44 Boo, QZ Car, SZ Cam, 
GT Mus) there were also the analysis of V1061 Cyg (combining the light c urve analysis, 



radial velocity analysis, light-time effect and Hipparcos measurements, see iTorres et al. 



20061). papers wher e radi al velocity measurements and astrometry were comb ined (see 



Muterspaugh et al.l (l2006l ) for the solution of the LITE system V819 Her, or iGudehus 
fl200lh "or IX Cas), the paper on HIP 50796 combining the radial- velocity measurements 
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with the Hipparcos abscissa data (see iTorred l2006l ) . or the paper on 6 Lib comparing the 



results from the period anal ysis, hght-curve analys is, spectral analysis, radio emission and 



astrometry, respectively; see iBudding et al.l (120051 ) 



Such a combined analysis is very important, and the individual methods could be 
tested. Their independent results have to be in agreement with each other. The method 
presented here is also the combination of the two independent methods into the one joint 
solution. It was never been done before in this way. Various modifications were presented, 
but most similar was the analy sis of the LITE tog ether with the Hipparcos observations. 



and the absolute astrometry by iRibas et al.l (120021 ) . 



The code itself is presented in section 12.91 It could be downloaded from the web sites 
and it is ready to be used. Short description of the code is presented and also the brief 
manual is available. Only slight modifications of the algorithm are necessary before the 
first run of the code. The numerics and the computing time required for the code strongly 
depends on the initial parameters and the input data (their quality and the size of the 
data set), but it could be slightly improved, see section [231 

A few eclipsing binaries were studied in this thesis. Detailed analysis was performed for 
QS Aql, VW Cep, ( Phe, V505 Sgr, HT Vir and V2388 Oph. These systems have relatively 
best coverage both in astrometry and LITE variation. This is the crucial part, as one can 
see in the case of VW Cep, which is the most suitable system for the simultaneous analysis. 

If precise measurements and good coverage of at least one period of the distant body 
in both methods are available, the presented method is very powerful and the parameters 
of the distant-body orbit could be derived very precisely. Even the distance of the system 
could be computed with high confidence level. The limiting factor is mainly the coverage 
of the orbit. In most of the cases the orbit was not covered sufficiently with data. Also 
almost all the systems included in the catalogue in chapter [4. 71 have only poor coverage of 
the orbit by data in both methods. 

VW Cep: The case where both methods have relatively best coverage of the orbit 
is VW Cep. In this case the resultant parameters of the third body satisfies the limit 
for the luminosity, and also the systemic velocity variations coincide with our hypothe- 
sis. New results are comparable with the previous ones. An additional fourth body was 
introduced to describe the long-term variation in times of minima, as well as in radial 
velocities. The system is probably more complicated than was assumed (chromospheric 
activity cycles, stellar spots and flares), and it was decided to explain only the most pro- 
nounced effects in the O — C diagram. Using the combined approach it is possible to 
derive the parallax to VW Cep more precisely than in any other previous papers, result- 
ing in TT = (35.85 ± 0.37) mas. The two different approaches (LITE3 + LITE4 and 
LITE3 -|- mass transfer) were used and their results compared. Both approaches lead to 
approximately the same results both in astrometry and times-of-minima analysis. The 
simultaneous analysis is able to describe the sy stem in its complexity and one has to dis- 
agree with the result by PribuUa et al. ( 20001 ). that the astrometric orbit could not be 



identified with the LITE3 variation from the O — C diagram. As one can see, our new 
results are in agreement with each other without any problems. On the other hand only 
further observations of this system will decide which approach (Solution I. or II.) is the 
right one. Two new times of minima were observed at Ondfejov Observatory. 



QS Aql: In the case of QS Aql the parameters of the distant-body orbit were mainly 
derived from the LITE analysis, because the coverage of the astrometric orbit is very poor 
and the old data are not very reliable. Due to this difficulty, the inclination of the orbit 
could not be derived precisely. The computed value of the inclination is quite low and 
the error quite high. Low inclination dictates high mass of the third body (but with large 
errors). The new derived mass of the third body is in contradiction with the previous 
photometric analysis, but one can get a consistent result within the error of this value. 

C Phe: The system ( Phe displays an apsidal motion together with the LITE and this 
explanation fits the O — C residuals quite well. This is the first time when the apsidal 
motion together with the LITE hypothesis were applied to this system. The astrometric 
analysis of ( Phe is complicated due to the fact that the period of the third body orbit is 
circa 3 times longer than the interval covered by the data. The time span of the minima 
measurements is even worse, only about one fourth of the orbit is covered. On the other 
hand the powerful combined analysis was able to estimate all of the parameters of the 
third-body orbit precisely. This approach lead to the period of the third body about 
220 years and the parameters of such a body, its predicted mass and the spectral type is 
in an excellent agreement with the previous photometric analysis. 

V505 Sgr: V505 Sgr is the system, where both the astrometry and also recent times 
of minima observations deviate from the predicted trend. Despite the fact the third body 
was detected more than 20 years ago, it could be even observable in the spectrum of 
the system, the complex figure of the system is still missing. The new result from the 
combined approach is in contradiction with the previous results from photometry and also 
spectroscopy. It indicates that the third body observable in spectra and light curve is 
different from the fourth body observable astrometrically. Only further detailed analysis 
would prove this hypothesis. 

HT Vir: The eclipsing system HT Vir is the case where the new value of mass of the 
distant body is about 2 times larger than one would expect. The distant component in the 
system is also a double and from the spectroscopy one is able to derive an upper limit for 
its mass. Regrettably, our new result is in contradiction with such a mass. This could be 
due to only a few times of minima observed in the linear part of the O — C diagram, new 
minima are needed in the next decades. Four new times of minimum light were observed. 
On the other hand the astrometric orbit is well-defined and almost whole orbit is covered. 

V2388 Oph: The last system is V2388 Oph, where two different approaches were used. 
The astrometric variation is rapid and since its discovery the third component has revolved 
a few times around the primary. On the other hand, the times of minima were obtained 
only rarely during the last decade. Due to this reason, the rapid change in order of 9 years 
is hardly detectable from the O — C diagram analysis and one could speculate about the 
inclination of the orbit. If the orbit's inclination is close to 180°, there could be no LITE 
evident in the O — C diagram, which was presented as another explanation and the mass 
transfer was suggested as an alternative explanation for the O — C diagram. Only further 
times of minima would prove or refuse this hypothesis. Another possible explanation is 
that the system is quadruple and the third body observable interferometrically is not the 
one which causes LITE. 

The final result is that the method itself is potentially very powerful but it is also very 



sensitive to the quality of the input data, especially if the method is used for determining 
the distance of these binaries. It can only be applied successfully in those cases where 
the astrometric orbit and the LITE in the O — C diagram are well defined by existing 
observations and lead to the approximately same parameters of the distant-body orbit. 
This is necessary condition, as one can see for example from the case V505 Sgr. 

The catalogue: The catalogue of other suggested systems for the prospective simulta- 
neous analysis with the introduced algorithm was presented in the chapter 14.71 The main 
purpose of the catalogue was to critically consider the potential objects for such a com- 
bined approach and from the eclipsing binaries in the spatially resolvable systems identify 
those, which are suitable for the introduced method. 

During the inspection of such systems there were found a few binaries which were often 
presented as eclipsing binaries, but which are in fact not. These are for example 18 Aql 
(=Y Aql) which is an ellipsoidal variable, o And which is photometrically variable, but it 
is not due to the eclipses, or V640 Cas which is probably also not an eclipsing variable. 
Also one misidentification of the EB was presented (HD 174932). 

On the other hand there were found a few systems which are the most suitable ones for 
the method presented here. Such systems are for example V348 And, V592 Per, V635 Mon, 
DN UMa, V831 Cen, ET Boo, or i Boo. 

Additional material is also available via the web page^. The code for computing the 
combined analysis could be downloaded together with the brief manual and instructions 
for the user. On the same web pages there are also the complete data files, which were 
used as the input files for the analysis. 



If 



http://sirrah.troja.mfF.cuni.cz/ ~ zasche/ 
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